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1.0 INTRODUCTION AND SUMMARY 


This user’s manual for the TETRA 2 program provides the requisite inputs and 
description to perform a steady state response analysis of a turbine engine. In 
addition, since TETRA 2 was developed with the dual capability of steady state 
response and the transient analysis of the original program, inputs for both 
options have been consolidated. This includes the common modeling procedures 
(and basic inputs) and the separate inputs unique to each option. Hence TETRA 2 
is a single program wherein the analytical and computational distinctions of the 
two options are transparent to the user. 

The transient option has all the capabilities of the original, including the 
later modules: 

(1) Flexible Bladed Disk 

(2) Squeeze Film Bearing 

In addition, the scope of the transient analysis has been increased so that it 
now contains the following enhancements initially developed for the steady state 
analysis: improved rub element (with the cubic nonlinear factor), structural 

damping capability (applies for physical connecting element types 1, 2, 4, and 
5), and new printout options. 

Also, changes in the coding of TETRA 2 were made to improve the efficiency of 
transient analyses. For most transient analysis runs, TETRA 2 is more than 
three times as fast as the earlier TETRA program. 

With the exception of the squeeze film bearing module, the added elements are 
also available in the steady state analysis. This exception was made because 
the squeeze film bearing is valid only for circular motion, whereas the steady 
state analysis is generally for elliptical orbits - which can exhibit large 
eccentricity or ellipticity. 

The driving motivation for the development of the steady state analysis was the 
need to calculate the steady unbalance response as a function of speed more 
directly. At the same time, it was needed to overcome the axi-symmetric 
limitations of current critical speed computer programs, and provide a 
computational method to calculate the realistic motions of actual engines. The 
non-isotropy in turbine support structures, such as those inherent in engine 
pylons and typical engine support hard points, require that the permissible 
engine motions are also non-isotropic. In other words, engine orbital motions 
are really elliptical, not circular. 

Current methods have required two analyses of critical speeds or steady state 
response - which are both axi-symmetric. In one, the structural stiffnesses 
correspond to the vertical plane stiffness; in the other, the horizontal plane 
stiffness. Additionally, this procedure does not truly represent the cross-axis 
(cross-planar) coupling intrinsic to gyroscopic moments on rotating and 
vibrating disks and rotors. 

TETRA 2 therefore was developed to overcome the shortcomings of standard steady 


1 


state analysis programs, at the same time providing new capabilities. These of 
course include the general nonlinear rub element (incorporating a deadband and a 
hardening or softening spring), the flexible bladed disks, and the modular 
design of the computer program. The latter allows the user to develop other 
connecting element modules or extend the basic structural components to 
accommodate more than four structural elements. 

The consolidation of the transient and the steady state options provide one 
computational method to perform most dynamic analyses of realistic turbine 
engine structural systems. Even though emphasis is on turbine engines, TETRA 2 
can also be used to perform dynamic analysis of non-turbine structures such as 
aircraft, buildings, or other structures which can be modeled with component 
structural components (and normal modes) connected by connecting elements. 

The charts presented in figures 1-1 and 1-2 graphically summarize the 
capabilities and operations of TETRA 2. The latter flow chart delineates TETRA 
2 in its two modes or options. The basic structural modeling and input 
information are common to both, including those for the connecting elements and 
other modules. Branching to either option is governed by the natural differences 
in the analyses: the transient option input/output include those that depend on 
time, such as rotor speed history, blade loss onset, number of time intervals, 
etc.; the steady state analysis input/output are in the frequency domain so that 
this branch requires information dependent on excitation frequency such as: 
frequency sweep, initial iteration solution trials, and others. 

Finally, due to the voluminous output of the transient solution, and in the case 
of critical speed calculations, of the steady state, a plot file output is 
recommended. This avoids the extensive printed output, at the same time 
providing quick and easily interpretable presentation of the results. Such a 
plot file is especially indispensable for other types of output processing, such 
as animation or FFT analysis of transient data, and the frequency response plot 
presentation of the steady state output. 

This volume contains the updated inputs for TETRA 2, combining the transient and 
steady state options. In addition explanations of the input variables and 
output are also included. 
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FIGURE 1-2 TETRA 2 SCHEMATIC DIAGRAM 
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2.0 TETRA 2 INPUT 


2.1 Input Restart File 

Transient analysis TETRA runs can be restarted using a restart file from an 
earlier TETRA run. A transient analysis run starts at time = 0 and marches 
forward in time. To begin a transient analysis run at a restart time (greater 
than 0), generalized coordinates and other data is needed for the restart time, 
and the generalized coordinates are also needed for the two time steps prior to 
the restart time. The restart file is the means of communicating this data 
calculated by the earlier TETRA run to the restarted TETRA run. For more 
information on restarting TETRA transient analysis runs, see the discussion of 
the type 1-1 input sheet. 

Steady state analysis runs cannot be restarted using a restart file from an 
earlier TETRA run. Restarting is unnecessary for a steady state analysis run, 
because steady state analyses can be started at any forcing frequency, and no 
knowledge of previous forcing frequencies that may or may not have been run is 
needed. The only possible exception to this is that the user might want to input 
an initial guess for the rub element relative displacements (the unknowns in the 
iterating equations) which was obtained from the printout at a given forcing 
frequency of a previous TETRA run. However, if this is the case the user can 
easily enter this initial guess in the primary input file (input sheet 1-3). 
Furthermore, an initial guess is not required. For more information regarding 
an initial guess, see the discussion of the type 1-3 input sheet. 

The input restart file is read in using file code 22. Therefore, if the run is 
a restart run, the user must include a line or two in his or her JCL (job 
control language) identifying the restart file and associating it with file code 


The restart file is an unformatted sequential file. It is a short file and 
takes up very little disk space. 
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2.2 Primary input file 

The primary TETRA input file consists of namelist input. The file code for the 
primary input file is the standard 05 file. See section 2.3 for a discussion of 
the namelist input sheets for this file. 

As in the original TETRA program, the namelist input is divided into sections of 
which there are four types named LIST1, LIST2, LIST3, and LIST4 (must be 
arranged in that order in the input file). Each section must begin with a $ 
followed immediately by the section name (that is, $LIST1, SLIST2, $LIST3, or 
SLIST4). Following this comes the input variable names and their values (which 
may be entered in any order). On the Honeywell computer, each section must end 
with a $. On the Cray computer, either a $ or a SEND can be used to signal the 
end of a section. The namelist input should begin in column 2 of each line. 

The LIST1 namelist input is for the overall input and input for points not 
located on the modal subsystems. There must be one and only one LIST1 section. 
Next come the LIST2 namelist sections. One LIST2 section is required for each 
modal subsystem (there must be at least one modal subsystem and at most 13). 

The LIST3 namelist sections (if any) come next. There is one LIST3 section for 
each physical connecting element. There may be zero physical connecting 
elements. The current upper limits are a maximum of 5 type 1 physical 
connecting elements, 10 of type 2, 10 of type 3, 4 of type 4, 25 of type 5, and 
5 of type 6. (Incidentally, the upper limits for the program variables may 
change in the future as the program evolves. A listing of the current upper 
limits is included in the printed output from a TETRA run just following the 
listing of the input file, except when the short form of the printed output is 

requested by entering IOUT * 0 on the type A input sheet). Finally, there must 

be one and only one LIST4 section. The LIST4 section includes restart input, 
time integration input, print and plot multiples, steady state nonlinear 
analysis input, rotor speed and rate input, load input, and plot file input. 

For a sample input file listing, see section 5.3.1. A listing of the input file 
for the demonstrator model is imbedded near the beginning of the printed output 

for the demonstrator model in this section. 
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2.3 TETRA 2 Input Sheets 


The TETRA 2 program has 44 namelist input sheets. However, usually many 
of the input sheets are not needed, and the user just includes those input 
sheets that apply for his or her model. Unless otherwise stated on the input 
sheet, the input data is common to both transient and steady state analyses. If 
the input is specific to a transient or to a steady state analysis, this is so 
noted on the input sheet. 

Examining the TETRA 2 input sheets, one can see that about 80% of the TETRA 
input (for instance, the usually lengthy modal subsystem input) is common to 
both transient and steady state analysis runs. Another 10% of the TETRA 2 input 
pertains specifically to a transient analysis run. Most of the input in these 
two categories was present in the earlier TETRA program (see references 1 and 
2). Finally, about 10% of the TETRA input is specific to a steady state 
analysis run. This input is all new with TETRA 2, of course, since the earlier 
TETRA program could only perform transient and not steady state analyses. 

Three enhancements were added to the TETRA 2 input that were made to apply for 
either transient or steady state analyses. These are: the rub element cubic 

nonlinear factor AK (on type F input sheet), structural damping frequency option 
ISF (on type A input sheet), and printout option IOUT (on type A input sheet). 

The TETRA 2 input was set up such that so that the same input file could easily 
be used for either a transient or steady state analysis. Simply include in this 
input file both the input that is specific to a transient analysis and the input 
that is specific to a steady state analysis. When a transient analysis is run, 
TETRA 2 just ignores the input that is specific to a steady state analysis, and 
likewise, when a steady state analysis is run, TETRA 2 ignores the input that is 
specific to a transient analysis. Then you need only change the value of input 
variable NTYPE (type A input sheet) to tell TETRA 2 whether the run is for a 
transient or steady state analysis. A transient analysis is the default if 
input variable NTYPE is not included. 

It is emphasized that the input for TETRA 2 is entirely consistent with that of 
the earlier TETRA program. In fact, if a user has an input file that he or she 
ran on the original TETRA program, that input file, without any modifications, 
will still work fine for TETRA 2, producing the same transient analysis run 
results as were produced by the original TETRA program. If the user then wants 
to use this input file for a steady state analysis run, he or she will have to 
add the variable NTYPE » 1 (type A input sheet) to tell the program that the run 
is for a steady state analysis and a small amount of additional input that is 
needed for a steady state analysis. 

Many of the input sheets to the earlier TETRA program remain the same for TETRA 
2. However, several input sheets to the earlier TETRA program were modified, 
and several others are new to TETRA 2. Although lots of this material is 
redundant of the earlier TETRA manuals (references 1 and 2), it was decided to 
discuss all the TETRA 2 input sheets in this manual, so that the user would only 
have to refer to one manual regarding the input sheets. 
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2.3.1 Overall Input (Type A Input Sheet) 

Except for the possible presence of several lines of JCL (job control language), 
the first line on the input file must be $LIST as shown on the type A input 
sheet. This line signals the beginning of the LIST1 namelist input section, 
which is comprised of type A and type B input sheets. Each input file must have 
one LIST1 namelist section. 

Next, the variables NAME, ADDRESS, IDENT1, and IDENT2 are used to enter the 
user’s name, address (building, mail drop, bin number, etc.), and two lines of 
description for the run. Up to 60 characters may be included between the 
apostrophes for each of these variables. These variables are not required, so 
one or more of them may be omitted if desired. If present, the information 
entered via these variables is reproduced on the first page (title page) of the 
printed output. 

Following this, the analysis type flag is entered. If NTYPE » 0, or if this 
variable is omitted, a transient analysis is performed. If NTYPE = 1, a steady 
state analysis is performed. 

Next comes the structural damping frequency option ISF. Structural damping is an 
enhancement which was added to TETRA 2 and was not present in earlier versions 
of TETRA. This feature can be used either for transient or steady state 
analyses. Structural damping may be used for any of the linear (type 1, type 2, 
type 4, or type 5) physical connecting elements, but is not available for the 
nonlinear physical connecting elements (type 3 or type 6). 

The presence or absence of the ISF input variable is not what determines whether 
structural damping is present, so if the user does not want structural damping, 
it does not matter whether this variable is entered or not or what value it is 
given. Rather, the user signals that he or she wants structural damping for a 
given physical connecting element by entering a non-zero value for QELEM 
(Q-factor) and by omitting QFREQ (the selected frequency) on the input sheets 
for that element (see the input sheets for the physical connecting elements, 
types D-l through H-3). 

If the user has so indicated that structural damping is desired for a particular 
element, the program then uses the entered Q-factor and the frequency defined by 
the ISF value to calculate the damping. If ISF equals 1 or is omitted, then the 
independent rotor speed is used regardless of whether the run is a transient or 
a steady state analysis. (If you don’t know what is meant by the independent 
rotor, read the discussion of the type J-l and J-2 input sheets). However, if 
the user enters ISF - 2, the independent rotor speed is still used if the run is 
a transient analysis run, but the steady state forcing frequency is used if the 
run is a steady state analysis run. (For steady state runs, the steady state 
forcing frequency might be the same as the independent rotor speed, but is 
different from the independent rotor speed when defined by the BW, FW, and SW 
variables on input sheet L-2). 
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Page of 

NAMELIST 
Type A 


IEJBA 

OVERALL INPUT 


Enter the following name, address, and run identification input (maximum of 60 
characters between apostrophes): 

SLIST1 

NAME-’ 

ADDRESS * , 

IDENTIC" * , 

IDENT2=’ ’ , 


Enter the analysis type flag NTYPE: 

NTYPE= , 0 - Transient analysis (assumed value) 

1 = Steady state analysis 


Enter the structural damping frequency option ISF, which defines the frequency 
to be used along with the Q-factor to calculate the structural damping (if 
present) for type 1, type 2, type 4, or type 5 physical connecting elements. 
ISF applies for each physical connecting element for which the user inputs a 
non-zero value for QELEM (Q-factor) and omits QFREQ (the selected frequency) 
(see the input sheets for the physical connecting elements). 

ISF- , 1 ■ Independent rotor speed used for both steady state and transient 

analyses (assumed value) 

2 * Steady state forcing frequency used for steady state analyses and 
independent rotor speed used for transient analyses 


Enter the printout option IOUT, which determines what output is printed. 

I0UT= , 0 ■ Short form of the printout 

1 - Standard form of the printout 

2 - Long form of the printout 
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If the user inputs both QELEM (the Q-factor) and QFREQ (the selected frequency) 
on the input sheets for the given physical connecting element, then the program 
uses the entered Q-factor and entered frequency to calculate the damping for the 
element (overriding the frequency defined by the ISF variable). If this is the 
case, the damping is non-structural , since the constant frequency defined by 
QFREQ is used regardless of variations in the independent rotor speed or the 
steady state forcing frequency. The ability to enter non-structural damping in 
this manner was present in the earlier versions of TETRA. 

The final variable entered on the type A input sheet is the printout option 
IOUT. Three possible values (0, 1, or 2) may be entered for this variable, but 
the value of 1 (standard form of the printout) usually is best (1 is also the 
default if this variable is omitted). A value of 1 provides a formatted listing 
of the input (very useful in debugging one’s input) and provides all the 
information that the user is normally interested in for the requested time steps 
(for a transient analysis) or the requested forcing frequencies (for a steady 
state analysis). This includes the physical displacements, physical forces, and 
all other essential data. 

If the user enters IOUT = 0, the short form of the printout is used. The short 
form of the printed output contains what the standard form of the output does 
minus the verbatim listing of the input file, the listing of the current program 
limits for various variables, and the formatted listing of the input as 
interpreted by the program. With the short form of the output, the title page 
is printed, followed by the same output for each time step (if a transient 
analysis) or for each forcing frequency (if a steady state analysis) that is 
printed for the standard output form. 

If the TETRA program detects any input errors, it prints an error message 
exactly at the point in the formatted listing of the program input that caused 
the error. This is highly useful to the user in figuring out what was wrong. 

If the user inputs IOUT » 0, the error message is still printed out, but since 
the formatted listing of the program input has been suppressed, the user 
sometimes is left scratching his head trying to figure out exactly where in the 
input the error occurred. Thus, it is recommended that the user not enter IOUT 
= 0 until after the program has been successfully run and the input thoroughly 
verified. Once the user is sure that the input is okay, setting IOUT » 0 will 
decrease the number of pages of output appreciably. 

If the user enters IOUT * 2, the long form of the printout is used. The long 
form of the output contains everything that the standard form of the output 
contains plus additional information. Usually, however, the additional items 
printed are of no interest to the user. Also, the long form of the output adds 
appreciably to the number of pages of printed output, especially for a steady 
state analysis run. 

The following items are added to the standard printout for transient analysis 
runs when the long form of the output is requested: more detailed 

information about the flexible bladed disk(s) (if present) (including 
printout of the mass matrix for the rotor that includes the flexible bladed 
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disks and for the flexible bladed disks), the (mode, point, direction) trios that 
contribute to the generalized forces, the generalized forces, generalized 
velocities, generalized displacements, generalized accelerations, and a listing 
of at least part of the output plot file (lists the first 11 time steps on the 
plot file for a transient analysis run). The listing of the output plot file 
would be helpful to someone writing a postprocessor program to plot the TETRA 
results. 

The following items are added to the standard printout for steady state analysis 
runs when the long form of the output is requested: more detailed 

% information about the flexible bladed disk(s) (if present) (including the terms 

of the global mass and global velocity matrices due to the flexible bladed 
disks), the (mode, point, direction) trios that contribute to the generalized 
forces, the difference in the mode shapes at the I and J ends for each rub (type 
3) element, a listing of at least part of the output plot file (lists the first 
11 forcing frequencies on the plot file for a steady state analysis run), 
several matrices that are calculated for each of the linear (type 1, 2, 4, or 5) 
physical connecting element and each gyroscopic element, the global matrices, 
various other matrices and vectors, the generalized forces, the generalized 
displacements, the parameters used in the iterating equations, the sine and 
cosine components of the physical displacements, the sine and cosine components 
of the physical velocities, the sine and cosine components of the modal forces, 
the sine and cosine components of the physical connecting element forces, and 
the sine and cosine components of the gyroscopic element forces. Note that 
while the standard form of the output omits the sine and cosine components of 
the physical displacements, physical velocities, modal forces, physical 
connecting element forces, and gyroscopic element forces, the standard form of 
the output does print out the magnitudes and phase angles for these quantities. 
The magnitudes and phase angles (not the sine cosine components) are also what 
gets written to the plot file (usually the user is just interested in plotting 
the magnitudes). 

Two other input variables also have significant effect on the length of the 
printed output. These are the print multiples (IPRMUL for a transient analysis 
and JPRMUL for a steady state analysis). These variables are entered on input 
sheet 1-1. See the discussion of the type 1-1 input sheet for details. 
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2.3.2 Physical Points Not Located On Modal Subsystems 
(Type B Input Sheet) 

This input sheet defines the boundary conditions of the problem. The points 
which the user desires to ground are defined here. If the code is specified as 
fixed (code = 1), the point is fixed in all 6 directions. If the code is 
specified as free (code = 0), the point is free to move all 6 directions. The 
points defined on this sheet may not be located on the modal subsystems but can 
only be located at ground or at the junctions of link and engine support 
elements. 
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TETRA 


Page of 

NAMELIST 
Type B 


PHYSICAL POINTS NOT LOCATED ON MODAL SUBSYSTEMS 


Z (vertical direction) 


y (horizontal direction) 


x (axial direction, positive forward) 


V 

PP= 

CODE 
0® FREE 
1 =FIXED 


* These points include ground points and points between link and engine 
support elements. It will be noted that points on modal subsystems 
cannot be attached directly to ground, but can be attached to ground 
through physical connecting elements. 


POINT NUMBER 
NOT ON MODAL 
SUBSYSTEM * 


COORDINATES RELATIVE 
TO GLOBAf 
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2.3.3 Modal Subsystem Input (Type C-l through C-16) 


Modal subsystem input is accomplished using the LIST2 namelist section. One 
LIST2 namelist section is required for each modal subsystem. The number of 
modal subsystems is counted by the program automatically. 

First, a discussion of the input that is common for all the modal subsystems. 

The user may input a title (this is optional) for each subsystem using the TITLE 
input variable. Up to 60 characters may be enclosed within the apostrophes for 
the title. 

The subsystem number (ISUB input variable) is required for all modal subsystems. 
The subsystem number must be an integer between 1 and 13. The subsystems are 
not arbitrary - each subsystem number represents a different type of modal 
subsystem with differing degrees of freedom (see table 2-1). The user chooses 
which of the modal subsystems to use, and can arrange these modal subsystems in 
any order. The user must have at least one modal subsystem. 

Note that subsystems 1 through 3 represent "rotor 1" and subsystems 4 through 6 
represent "rotor 2" (see table 2-1). Of course, the user’s model may have only 
one rotor, in which case either rotor 1 (subsystems 1 through 3) or rotor 2 
(subsystems 4 through 6) could be used. If the engine is a dual spool engine, 
however, the user might want to use both "rotor 1" and "rotor 2" in his model. 

In this instance, "rotor 1" might represent, say, the low pressure system (fan, 
low pressure compressor, shaft, and low pressure turbine), and "rotor 2" might 
represent the high pressure system (high pressure compressor, shaft, and high 
pressure turbine). Or vice-versa - the choice is up to the user. The user has a 
great deal of flexibility and with these subsystems can model virtually any jet 
engine configuration. 

Next, the input variables XREF, YREF, ZREF, and PTS are used for almost all of 
the subsystems. Specifically, these variables are used for the rotor, case, and 
pylon subsystems (subsystem numbers 1 through 11), but are not used for the 
flexible bladed disk subsystems (subsystems 12 and 13). If the subsystem 
reference point does not coincide with the global origin, the user must input 
the coordinates of the subsystem reference point (XREF, YREF, AND ZREF) input 
variables. Then, the user enters the points of the subsystems and their 
coordinates with respect to the subsystem reference point via the PTS input 
variable. Each subsystem must have at least one point. The current upper 
limits are 10 points for the rotor' subsystems (subsystem numbers 1 through 6), 

20 points for the case subsystems (subsystem numbers 7 through 10), and 10 
points for the pylon subsystem (subsystem number 11). 

The following sections detail the input that is specific to a particular type 
of subsystem. 


Table 2-1. Modal Subsystems. 


SUBSYSTEM 

NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 
13 


NUMBER OF 
DEGREES OF 


DESCRIPTION FREEDOM 

Rotor 1 Vertical Plane 2 

Rotor 1 Horizontal Plane 2 

Rotor 1 Rigid Body 5 

Rotor 2 Vertical Plane 2 

Rotor 2 Horizontal Plane 2 

Rotor 2 Rigid Body 5 

Case Vertical Plane 2 

Case Horizontal Plane 2 

Case Rigid Body 6 

Case Torsional 1 

Pyl on 3 

Flexible Bladed Disk Number 1 2 

Flexible Bladed Disk Number 2 2 
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2.3.3. 1 Vertical and Horizontal Plane Subsystems (Type C-l through C-3) 

Input sheets C-l through C-3 apply for the vertical and horizontal plane 
subsystems, that is, subsystem 1 (rotor 1 vertical plane), subsystem 2 (rotor 1 
horizontal plane), subsystem 4 (rotor 2 vertical plane), subsystem 5 (rotor 2 
horizontal plane), subsystem 7 (case vertical plane), and subsystem 8 (case 
horizontal plane). The first input variables for the vertical and horizontal 
plane subsystems are TITLE, IS.UB, XREF, YREF, ZREF, and PTS. These variables 
were discussed in section 2.3.3. 

Next, the frequency w, potential energy PE, Q-factor Qf, and mode type code are 
entered for each mode using the XMODES input variable. There must be at least 
one mode entered. The current upper limit is 15 modes for the rotor subsystems 
(subsystem numbers 1 through 6), 30 modes for the case subsystems (subsystem 
numbers 7 through 10), and 30 modes for the pylon subsystem (subsystem number 
11 ). 

The input frequency to and input potential energy PE are used to calculate the 
generalized mass ZM for the mode as follows: 


2PE 



The Q-factor Qf determines the modal damping for the mode (provided that the 
mode type code equals 0) and is defined as follows: 

1 

Qf - 

2(C/C c ) 

where C is the damping coefficient and C c is the critical damping coefficient. 

If the user wishes to neglect damping, Qf=0 should be inputted, which is a 
signal to the computer to neglect damping. 

The input mode type code must be 0 (flexible) or 1 (rigid body). The code 1 is 
used for "soft spring" rigid body modes, which are used to approximate the true 
rigid body modes. If this code equals 0 (flexible), the modal stiffness ZK and 
the modal damping ZC are calculated as follows: 

ZK 

ZK = 2PE ; ZC = 

wQf 

However, if the mode type code equals 1 (rigid body), then the modal stiffness 
and modal damping are both se+ to 0, regardless of what was input for PE and Qf. 

Finally, the mode shape input (translation, slope, shear, and moment) must be 
entered using the VH input variable. Mode shapes must be entered for every mode 
.that was included in the XMODES input. 
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TETRA 


PAGE OF 

NAMELIST 
TYPE C-l 


MODAL SUBSYSTEM INPUT FOR VERTICAL AND 
HORIZONTAL PLANE SUBSYSTEMS 


THIS INPUT SHEET APPLIES FOR THE FOLLOWING MODAL SUBSYSTEMS: 


NUMBER 1 
NUMBER 2 
NUMBER 4 
NUMBER 5 
NUMBER 7 
NUMBER 8 


(ROTOR 1 VERTICAL PLANE) 
(ROTOR 1 HORIZONTAL PLANE) 
(ROTOR 2 VERTICAL PLANE) 
(ROTOR 2 HORIZONTAL PLANE) 
(CASE VERTICAL PLANE) 

(CASE HORIZONTAL PLANE) 


Z 



Z (VERTICAL DIRECTION) 



POSITIVE FORWARD) 


NOTE - THE PHYSICAL MODEL (MODAL SOURCE) FOR THE SUBSYSTEM IS ALWAYS A 


VERTICAL PLANE MODEL. 


T 

$LIST2 

TITLE* ’ 

ISUB= , SUBSYSTEM NUMBER (1,2, 4, 5, 7, or 8) 

ENTER THE COORDINATES (INCHES) OF THE SUBSYSTEM REFERENCE POINT: 
XREF= ,X COORDINATE (0 ASSUMED) 

YREF= ,Y COORDINATE (0 ASSUMED) 

ZREF= ,Z COORDINATE (0 ASSUMED) 


Maximum of 60 characters enclosed within apostrophes for the title. 
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LOCAL POINT NUMBER 


PAGE OF 

NAMELIST 

TETRA TYPE C-2 


MODAL SUBSYSTEM INPUT FOR VERTICAL AND 
HORIZONTM. PLANE SUBSYSTEMS (CONTINUED) 


ENTER THE COORDINATES OF THE TETRA POINTS. MAXIMUM 0. 10 
POINTS FOR SUBSYSTEMS 1,2, 4, AND 5. MAXIMUM OF 20 POINTS FOR 
SUBSYSTEMS 7 AND 8. Y and Z COORDINATES NORMALLY 0. 
n» NUMBER OF POINTS IN THE SUBSYSTEM. 


TETRA 

POINT 

COORDINATES (INCHES) 
SUBSYSTEM REFERENCE 

WITH RESPECT TO 
POINT 

NUMBER 


X Y 

Z 

w 

PTS = 

9 

9 

9 9 



9 

9 9 


9 

9 

9 9 


n , * • » 


ENTER THE VALUES BELOW FOR EACH MODE. MAXIMUM OF 
30 MODES. 

N=NUMBER OF SUBSYSTEM MODES. 


FREQUENCY 

POTENTIAL 

Q-FACTOR 

MODE TYPE 
0*FLEXIBLE 

(RPM) 

ENERGY 


1-RIGID BODY 




XMODES « 
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TETRA 


PAGE OF 

NAMELIST 
TYPE C-3 


MODAL SUBSYSTEM INPUT FOR VERTICAL AND 
HORIZONTAL PLANE SUBSYSTEMS (CONTINUED) 


ENTER THE REQUIRED MODE SHAPES BELOW. 
n-NUMBER OF POINTS IN THE SUBSYSTEM. 
N=NUMBER OF SUBSYSTEM MODES. 


TRANSLATION T 

SLOPE e 

SHEAR V 

MOMENT M 

INCHES 

RADIANS 

POUNDS 

IN-LB 


VH(1,1,1)= 


C5 

1 

9 

9 

9 

9 

W 

2 

9 

9 

r 

9 

3 

9 

9 

9 

9 

H g 






sg 

• 

• 




g2 

• 

• 




O J 

J 3 

• 

« 

iL 

• 

• 

• 

s 

9 

9 

9 



1 

VH(1,1,2)= 

9 

9 

9 

9 

63 

2 


9 

9 

9 

9 

<N § 

3 


9 

9 

9 

9 

LOCAL MODE 
LOCAL POINT Nl 

• 

• 

• 

n 

• 

9 

9 

9 

9 


VH(1, 1,N)= 


cs 

1 

9 

9 

9 

9 

W 

* I 

2 

, 

9 

9 

9 

3 

9 

9 

9 

9 

H as 






as 

• 

• 




< 2 
8 j 

• 

• 

m 

• 




S 

• 

• 

• 





n 

* 

9 

9 

9 
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2. 3. 3. 2 Rigid Body Subsystems (Type C-4 through C-6) 

Input sheets C-4 through C-6 apply for the rigid body subsystems, that is, 
subsystem 3 (rotor 1 rigid body), subsystem 6 (rotor 2 rigid body), and 
subsystem 9 (case rigid body). The first input variables for the rigid body 
subsystems are TITLE, ISUB, XREF, YREF, ZREF, and PTS. These input variables 
were discussed in section 2.3.3. 

Next, input variables XC6, YCG, and ZCG are used to input the coordinates of the 
subsystem center of gravity with respect to the subsystem reference point (0 is 
the assumed value for these quantities if not input). 

In addition, the physical weight properties (WX, WY, WZ, WMIX, WMIY, and WMIZ) 
should be input, but only for those directions for which rigid body motion is to 
be considered. Generalized coordinates are assigned for each direction for 
which the physical weight property is non-zero, but are not assigned if the 
physical weight property was set to zero or omitted. Since the e x direction is 
always neglected for subsystems 3 and 6, WMIX should be set to zero or omitted 
for subsystems 3 and 6. At least one of the physical weight properties must be 
non-zero. 

Note that there are two ways of accounting for the rigid body modes, that is, 

by entering the "soft spring" rigid body modes with the vertical and horizontal 

plane subsystem input (setting the mode type code to 1) (input sheets C-l 
through C-3), or entering the rigid body modes with the rigid body subsystem 
input (input sheets C-4 through C-6). The user is cautioned not to enter these 
modes both places to avoid doubling the effect of these modes. Both ways of 

modeling the rigid body motion give equivalent results, as confirmed by test 

cases. However, it has been found that TETRA runs faster (and thus saves money) 
if the rigid body modes are included in the vertical and horizontal plane 
subsystems rather than the rigid body subsystems. Thus, it is recommended that 
the user use the vertical and horizontal plane subsystems in lieu of the 
rigid body subsystems for the rigid body modes. An exception to this guideline 
is if the user needs to account for rigid body motion in the X (axial) or 5 X 
(torsional) direction, in which case the rigid body subsystem would be needed 
since the vertical and horizontal plane subsystems do not consider these 
directions. 
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TETRA 


PAGE OF 

NAMELIST 
TYPE C-4 


MODAL SUBSYSTEM INPUT FOR RIGID BODY SUBSYSTEMS 


THIS INPUT SHEET APPLIES FOR MODAL SUBSYSTEM 3 (ROTOR 1 RIGID BODY), MODAL 
SUBSYSTEM 6 (ROTOR 2 RIGID BODY), AND MODAL SUBSYSTEM 9 (CASE RIGID BODY). 


Z 



Z (VERTICAL DIRECTION) 



X (AXIAL 
DIRECTION) 


NOTE: MOTION IN THE 3X DIRECTION IS NOT CONSIDERED FOR THE ROTOR SUBSYSTEMS 

(3 AND 6). MOTION IN THIS DIRECTION CAN BE CONSIDERED FOR THE CASE SUBSYSTEM 
(9), HOWEVER. 

W 

$ 

$LIST2 

TITLE = , 

ISUB= , SUBSYSTEM NUMBER (3, 6, OR 9) 

ENTER THE COORDINATES (INCHES) OF THE SUBSYSTEM REFERENCE POINT: 

XREF= , X COORDINATE (O ASSUMED) 

YREF= , Y COORDINATE (O ASSUMED) 

ZREF= , Z COORDINATE (O ASSUMED) 


Maximum of 60 characters enclosed within apostrophes for the title. 


* 
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LOCAL POINT NUMBER 


PAGE OF 

TETRA NAMELIST 

TYPE C-5 

MODAL SUBSYSTEM INPUT FOR RIGID BODY SUBSYSTEMS (CONTINUED) 

ENTER COORDINATES OF TETRA POINTS. MAXIMUM OF 10 POINTS FOR 
SUBSYSTEMS 3 AND 6 AND MAXIMUM OF 20 POINTS FOR SUBSYSTEM 9. 

Y AND Z COORDINATES NORMALLY 0. 

QS NUMBER OF SUBSYSTEM POINTS. 


COORDINATES (INCHES) WITH RESPECT 
TO SUBSYSTEM REFERENCE POINT 



ENTER THE SUBSYSTEM CENTER OF GRAVITY COORDINATES 
RESPECT TO SUBSYSTEM REFERENCE POINT: 



X COORDINATE 
Y COORDINATE 
Z COORDINATE 


(0 ASSUMED) 
(0 ASSUMED) 
(0 ASSUMED) 


(INCHES) WITH 
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TETRA 


PAGE OF 

NAMELIST 
TYPE C-6 


MODAL SUBSYSTEM INPUT FOR RIGID BODY SUBSYSTEMS (CONTINUED) 


ENTER THE SUBSYSTEM PHYSICAL WEIGHT PROPERTIES FOR EACH OF THE SIX 
DIRECTIONS. IF MOTION IN A PARTICULAR DIRECTION IS TO BE NEGLECTED, 
ENTER ZERO FOR THE PHYSICAL WEIGHT PROPERTY IN THAT DIRECTION. 

(SINCE THE 0X DIRECTION IS ALWAYS NEGLECTED FOR SUBSYSTEMS 3 AND 6, 
WMIX SHOULD ALWAYS BE SET TO ZERO OR OMITTED FOR SUBSYSTEMS 3 AND 6) 
^ LEAST ONE OF THE PHYSICAL WEIGHTS MUST BE NON-ZERO. 




X 

wx= 

OSes 

B 

Y 

WY- 

M 

s 

w 

z 

wz= 

MOW 
W M CL. 

ex 

WMIX = 

> w o 
E ^ as 

ei 

M 

ey 

WMIY= 

CU 

o 

0Z 

WMIZ = 


WEIGHT, POUNDS (0 ASSUMED) 
WEIGHT, POUNDS (0 ASSUMED) 
WEIGHT, POUNDS (0 ASSUMED) 

POLAR MOMENT OF INERTIA, LB- IN 2 
MOMENT OF INERTIA ABOUT Y AXIS, 


(0 ASSUMED) 

LB- IN 2 (0 ASSUMED) 


MOMENT OF INERTIA ABOUT Z AXIS, LB-IN' 5 (0 ASSUMED) 
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2. 3. 3. 3 Case Torsional Subsystem (Type C-7 through C-9) 


Input sheets C-7 through C-9 are used for the case torsional subsystem, which 
subsystem 10. The first input variables for this subsystem are TITLE, ISUB, 
XREF, YREF, ZREF, and PTS. These input variables were discussed in section 
2.3.3. Next comes input variable XMODES, which was discussed in section 
2. 3. 3.1. 

The final input variable for the torsional subsystem is variable TOR, which is 
used to enter the mode shapes (twist and moment). Mode shapes must be entered 
for every mode that was included in the XMODES input. There must be at least 
one mode. The current upper limit for the case subsystems is 30 modes. 



TETRA 


PAGE OF 

NAMELIST 
TYPE C-7 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 10 
(CASE TORSIONAL SUBSYSTEM) 


Z (VERTICAL DIRECTION) 



V 

$ 

$LIST2 

TITLE= » 

ISUB=10, 

ENTER THE COORDINATES (INCHES) OF THE SUBSYSTEM REFERENCE POINT: 
XR£F= ,X COORDINATE (0 ASSUMED) 

YREF= ,Y COORDINATE (0 ASSUMED) 

ZREF- ,Z COORDINATE (0 ASSUMED) 


Maximum of 60 characters enclosed within apostrophes for the title. 
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TETRA 


PAGE OF 

NAMELIST 
TYPE C-8 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 10 
(CASE TORSIONAL SUBSYSTEM) (CONTINUED) 


ENTER THE COORDINATES OF THE TETRA POINTS. MAXIMUM OF 20 POINTS. 
X AND Y COORDINATES NORMALLY 0 
n*NUMBER OF POINTS IN THE SUBSYSTEM 


\ 

1 

PS 

w 

2 

it 

3 

i 

• 

6-« 

55 

• 

HH 

• 

2 

• 


• 

S 

• 

s 

• 

n 


TETRA 

POINT 

COORDINATES (INCHES) WITH RESPECT 
TO SUBSYSTEM REFERENCE POINT 

NUM3ER 

r x y i z 

V 


PTS = 

> * 

: ...» 

• : 



ENTER THE VALUES BELOW FOR EACH MODE. MAXIMUM OF 30 
N-NUMBER OF SUBSYSTEM MODES. 


MODES . 


FREQUENCY 

(RPM) 

POTENTIAL 

ENERGY 

Q-FACTOR 

MODE TYPE 
0=FLEXIBLE 
1“RIGID BODY 

V 

XMODES= 

t 

9 

9 

» 


9 

9 

9 

* 


9 

9 

9 

9 

• 

• 

• 

• 

• 

9 

9 

9 

9 


PS 

l 

W: 

2 

Si 

55 

3 

W 

• 

§ 

• 

s 

• 

iji 

• 

s 

• 

O 

<4 

N 
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TETRA 


PAGE OF 

NAMELIST 
TYPE 09 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 10 
(CASE TORSIONAL SUBSYSTEM) (CONTINUED) 


ENTER THE REQUIRED MODE SHAPES BELOW. 
n-NUMBER OF POINTS IN THE SUBSYSTEM. 
N-NUMBER OF SUBSYSTEM MODES. 



TWIST 9 

MOMENT 


RADIANS 

IN-LB 
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2. 3. 3. 4 Pvlon Subsystem (Type C-10 through C-12) 

Input sheets C-10 through C-12 are used for the pylon subsystem, which is 
subsystem 11. The first input variables for this subsystem are TITLE, ISUB, 
XREF, YREF, ZREF, and PTS. These input variables were discussed in section 
2.3.3. Next comes input variable XMODES, which was discussed in section 
2.3.3. 1. 

The final input variable for the pylon subsystem is variable PYL, which is used 
to enter the mode shapes (translation in the X, Y, and Z directions). Mode 
shapes must be entered for every mode that was included in the XMODES input. 
There must be at least one mode. The current upper limit for the pylon 
subsystem is 30 modes. 



TETRA 


PAGE OF 

NAMELIST 
TYPE C-IO 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 11 
PYLON SUBSYSTEM 


Z Z (VERTICAL DIRECTION) 



W 

$ 

$LIST2 

TITLE* ' 

ISUB=11, 

ENTER THE COORDINATES (INCHES) OF THE SUBSYSTEM REFERENCE POINT: 
XREF= , X COORDINATE (0 ASSUMED) 

YREF= , Y COORDINATE (0 ASSUMED) 

ZREF= , Z COORDINATE (0 ASSUMED) 


Maximum of 60 characters enclosed within apostrophes for the title. 
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PAGE 


OF 


TETRA 


NAMELIST 
TYPE C-ll 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 11 
(PYLON SUBSYSTEM) (CONTINUED) 

ENTER THE COORDINATES OF THE TETRA POINTS. MAXIMUM OF 10 POINTS. 
nsjJUMBER OF POINTS IN THE SUBSYSTEM 


TETRA 

POINT 

NUMBER 

COORDINATES (INCHES) WITH RESPECT 
TO SUBSYSTEM REFERENCE POINT 

X 

Y 

z 


W 

PTS = 


05 

1 

2 

9 

9 

’ 

i 

9 — 

} 

w 

3 

9 

9 

9 

-9 

E-» 

• 

• 

• 

• 




S3 

M 

• 

• 




2 


• 




iJ 

• 





§ 

• 

• 




O 

• 

n 

• 

9 


9 


ENTER THE VALUES BELOW FOR EACH MODE. MAXIMUM OF 30 MODES. 
N=NUMBER OF SUBSYSTEM MODES. 


FREQUENCY 

(RPM) 

POTENTIAL 

ENERGY 

Q-FACTOR 

MODE TYPE 
0=FLEXIBLE 
1=RIGID BODY 

XMODES = 

9 

- 1 

-9 

9 

9 

. J 

i 9 


9 

_ _J 

> 9 

-9 

• 

• 

• 

• 


7 

.9 


05 

W 

1 

jj§ 

2 


3 

W 

• 

P 


O 

• 

s 



# 

< 

• 

o 

o 

N 
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LOCAL MODE H LO CAL MODE 2 LOCAL MODE 1 

LOCAL POINT NUMBER LOCAL POINT NUMBER LOCAL POINT NUMBER 

-J..* • • to tO H p • ■ • • • OJ to H 0 CO to H 


TETRA 


PAGE OF 

NAMELIST 
TYPE C-12 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 11 
PYLON SUBSYSTEM (CONTINUED) ~ ' 


ENTER THE REQUIRED MODE SHAPES BELOW, 
as NUMBER OF POINTS IN THE SUBSYSTEM. 
N= NUMBER OF SUBSYSTEM MODES. 



PYL (1,1,2) = 



PYL (1,1, N) = 
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2. 3. 3. 5 Flexible Bladed Disk Subsystems (Type C-13 through C-16) 

Input sheets C-13 through C- 16 are used for the flexible bladed disk subsystems, 
that is, subsystem 12 (flexible bladed disk number 1) and subsystem 13 (flexible 
bladed disk number 2). For more details pertaining to the flexible bladed 
disks, see reference 2. 

A TETRA run can have 0, 1, or 2 flexible bladed disks. If two flexible bladed 
disks are present, they must be located on the same rotor. If the flexible 
bladed disks are located on rotor 1, then the rigid body modes for rotor 1 must 
be included in subsystems 1 and 2, and subsystem 3 must be omitted. If the 
flexible bladed disks are on rotor 2, then the rigid body modes for rotor 2 must 
be included in subsystems 4 and 5, and subsystem 6 must be omitted. Both the 
rotor vertical and horizontal plane subsystems should be included for the rotor 
which includes the flexible bladed disk(s) to account for coupling between the 
generalized coordinates. The physical weight and rigid moment of inertia 
properties of the flexible bladed disk(s) must be included in the applicable 
rotor vertical and horizontal plane subsystems. To account for gyroscopic 
effects, the polar moment of inertia of the flexible bladed disk must be 
entered on input sheet N at the flexible bladed disk center of gravity 
point. Rotor speed (input sheets J-l and J-2) must be included for the rotor on 
which the flexible bladed disk(s) are located. 

The first input variables are TITLE and ISUB. These input variables were 
discussed in section 2.3.3. The next input variable is ICG, the TETRA point 
number for the center of gravity of the flexible bladed disk. This point must 
be included in the applicable rotor vertical and horizontal plane input. 

Next, several flexible bladed disk properties are input, namely WTF (one 
diameter modal weight), QFAC (the modal Q-factor), WF (static frequency), XMU 
(modal tangential shear coefficient), SV (modal moment coefficient), and BETA 
(rpm-beta factor table). If damping is not desired, the modal Q-factor QFAC 
should be set to 0 or omitted, which is a signal to the program to neglect 
damping. The rpm-beta factor table must include at least 2 and a maximum of 10 
lines, with the table entries in order of increasing rpm. 

Following this comes the input for the local points on the flexible bladed disk. 
This includes input variables FBDPTS (the radius and polar angle for each local 
point), FBDMS (the static mode shapes, that is, tangential translation and axial 
translation -at the local points), and FBDS (the three modal stress components at 
the local points). Displacements and stresses are output for each of the local 
points. There may be zero local points if the user is not interested in 
flexible bladed disk displacements and stresses. The maximum number of local 
points permitted is 200 for each flexible bladed disk. 
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ORIGINAL TTOE » 
OF POOR QUALITY 


MODAL SUBSYSTEM INPUT FOR FLEXIBLE BLADED DISK SUBSYSTEMS 

Input aheeta C-13 through C-16 apply for Subayatea 12 (Flexible Bladed 
Diak Number 1) and Subayatea 13 (Flexible Bladed Diak Number 2). 

If both FBD Number 1 and FBD Nuaber 2 are preaent, they auat be located 
on the aaae rotor. If the FBD(a) are on Botor 1, then the rigid body 
aodea for Botor 1 auat be included in aodal Subayateaa 1 and 2, and 
aodal Subayatea 3 auat be omitted. If the FBD(a) are on Botor 2, then the 
rigid body aodea for Botor 2 auat be included in Subayateaa 4 and 5, and 
Subayatea 6 auat be omitted . Both the rotor vertical plane and horizontal 
plane aubayataaa ahould be included for the rotor which includea the 
FBD(a) to account for coupling between the generalized coordinatea. The 
phyaical weight and rigid aoaent of inertia propertiea of the FBD(a) auat 
be included in the applicable rotor vertical and horizontal plane aub- 
ayateaa. Botor apeed input (Input Sheets J-l and J-2) must be included 
for the rotor on which the FBD(s) are located. 

5 END 
SLIST2 
TITLE-' 

ISUB- __ 

ICG- 


Maximum of 60 characters enclosed within apostrophes for the title- 


I 

, SUBSYSTEM NUMBEB (12 or 13) 

,TETBA Point nuaber for the FBD center of gravity (Thia 
point auat be included in the applicable rotor vertical 
and horizontal plane aubayatea input) 


Page of 

NAMELIST 
Type C-13 
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Page of 

NAMELIST 
Type C-14 

MODAL SUBSYSTEM INPUT FOR FLEXIBLE BLADED DISK SUBSYSTEMS 

V 

WTF* , one-diameter modal weight WTf (lb) 

QFAC* , modal Q-factor (omit or set to 0 if no damping 

desired) 

WF- , static frequency Wf (hertz) 

XMU— , modal tangential shear coefficient. Mu (lb) 

SV— , modal moment coefficient, S v (in-lb) 

Note: To account for gyroscopic effects for the FBD, the polar moment of 

inertia of the FBD must be input on Sheet N at the FBD C.G. point. 

Input the following table for beta factors. Include at least 2 and a 
maximum of 10 lines. Entries must be in order of increasing rpm. 



B 

RPM 

FACTOR 


V 

beta- 
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__ page of 

ORIGINAL PAGE 15 NAMELIST - - 

OF POOR QUALITY Type C-15 


MODAL SUBSYSTEM INPUT FOR FLEXIBLE BLADE D DISK SUBSYSTEMS (Continued) 


Z (Vertical Direction) 



Enter the r and € coordinates for each of the subsystem points other than the 
center of gravity point. Maximum of 200 points for subsystem 12 and 200 
points for subsystem 13 other than the center of gravity points, n ■ number 
of subsystem points other than the center of gravity point. 




r > 0 

0* < * < 360* 


FBDPTS- 




Polar 

Radius 

Angle 

r 

C 

inches 

degrees 
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Page of 

NAMELIST 
Type C-16 


MODAL SUBSYSTEM INPUT FOR FLEXIBLE BLADED DISK SUBSYSTEMS (Concluded) 


For each of the aubaystem points other than the C.G. point enter the 
atatic (zero apeed) mode shapes (assumed to be 0 if not entered): 


FBD FBD 

Tangential Axial 

Translation Translation 

u (inches) v (inches) 


FBDMS - 



stress components (assumed to be 0 if not entered): 


01 

02 

03 

psi 

psi 

psi 




FBDS - 
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2.3.4. Physical Connecting Element Input (Type D-l through H-31 


TETRA 2 has capability for 6 types of physical connecting elements. These types 
are the general spring-damper element (type 1), the space link-damper element 
(type 2), the rub element (type 3), the engine support-links element (type 4), 
the uncoupled point spring-damper element (type 5), and the squeeze film damper 
element (type 6). All of these types may be used for either transient or steady 
state analysis runs with the exception of the squeeze film damper element, which 
may be used only for transient analysis runs. 

A new feature added to TETRA 2 is structural damping capability, which applies 
for the linear (type 1, 2, 4, and 5) physical connecting elements. This feature 
can be used either for transient or steady state analyses. As one can see from 
the input sheets, the input variables QELEM (q-factor) and QFREQ (selected 
frequency) may be used for any of the linear physical connecting elements. If 
both QELEM and QFREQ are entered, then damping is calculated based on the 
entered q-factor and entered selected frequency. This type of damping is 
non-structural , and was present in the earlier versions of TETRA. However, if 
QELEM is entered but not QFREQ, then damping is based on the entered Q-factor 
and either the independent rotor speed or the steady state forcing frequency, as 
determined by input variable ISF on the type A input sheet (see the discussion 
of the type A input sheet). This type of damping, which varies as the 
independent rotor speed or the steady state forcing frequency changes, is 
structural damping, and was not present in the earlier versions of TETRA. 
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2.3.4. 1 General Spring-Damper Element (Type D-l and D-2) 

This element is associated with two physical points located at arbitrary 
locations in global space. Each of these points is assigned six degrees of 
freedom, three translational displacements (or velocities), and three rotational 
displacements (or velocities). 

A full complement of stiffness and damping coefficients can be input to allow 
the modeling of fully coupled load paths. The units of the stiffness matrix 
coefficients are: Ib/in, lb/rad, in-lb/in, and in-lb/rad. The units of the 

damping matrix coefficients are Ib-sec/in, lb-sec/rad, in-lb-sec/in, and 
in-lb-sec/rad. 

Damping, if present, can be specified directly via damping matrix definition 
(non-structural damping), may be calculated from the input q-factor and selected 
frequency (non-structural damping), or may be calculated based on the input 
q-factor and either the independent rotor speed or the steady state forcing 
frequency as determined by input variable ISF on the type A input sheet 
(structural damping). See section 2.3.4 for further discussion of structural 
damping. 
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TYPE 1 PHYSICAL CONNECTING ELEMENT 



X (axial direction, positive forward) 


v 

$ 

$LIST3 
ITYPE=1 , 

ILEM= , element number 

I -end J-end 

point number point number 


STIFFNESS MATRIX DEFINITION 


I-END J-END 

GLOBAL DIRECTION DISPLACEMENT GLOBAL DIRECTION DISPLACEMENT 
~ I 2 H Tl |5 Ti 1 I 3 I 4 I 5 I 6 
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TYPE 1 PHYSICAL CONNECTING ELEMENT 
(GENERAL SPRING-DAMPER ELEMENT) (Continued) 


V 

I DAMP* 


1 “damping 
0*no damping 


(0 assumed) 


If I DAMP* 1 , complete the input for one of the following 
two options: 


Option 1: For damping based on Q-factor and selected 

frequency, enter the following (if QELEM is entered 
but not QFREQ, the frequency defined by the ISF 
q-factor variable on type A input sheet is used). 

QELEM* , 

frequency (hertz) 

QFREQ* , 


Option 2: For damping based on damping matrix definition, 

enter the following: 


DAMPING MATRIX DEFINITION 


I-END 

J-END i 

GLOBAL DIRECTION VELOCITY 

GLOBAL DIRECTION VELOCITY 


1 

2 

3 

EH 

5 

6 

n 

2 

3 

4 


6 

• 

Z 

°y 

m 

* 

9 z 

■ 

®x 

B 

6 y 

Y 

®z 

H 

®x 


G1 obal 

Dir. Forces 

DAMP = 

9 


9 

9 

9 

9 


1 

F, 


2 

F* 

9 

9 

9 

9 

9 

9 

O 

3 

F*/ 

9 

9 

> 

9 

9 

9 

Ui 

4 

F ez 

9 

9 

9 

9 

9 

9 

1 

5 

F ! 

hh 

— 

mm 

warn 

mm 

m 


6 

F 

■BG 

Bfi 

56 

56 

56 

■ 



0 X 








1 

HBB 

9 

9 



9 

9 


Q 

mm 


9 


9 


9 

o 

z 

MM 



9 


9 


9 

UJ 

1 

H 

Vsjbk 





9 


n 

H 









H 

ISBii 



9 





■■ 

HiH 

i 
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2. 3. 4. 2 Space Link-Damper Element (Type E) 


This element is used tc mcdel load paths which have no local moment capability. 
The load paths are associated with two physical points located at arbitrary 
locations in global space. Each of these points is assigned three translational 
degrees of freedom. The user inputs the cross sectional area and Young’s 
modulus, and TETRA 2 uses this information along with the coordinates of the two 
connecting points to calculate the stiffness matrix. 

Damping, if present, can be specified by the translational damping coefficient 
(non-structural damping), by the input q-factor and selected frequency 
(non-structural damping), and by the input q-factor and either the independent 
rotor speed or the steady state forcing frequency as determined by variable ISF 
on the type A input sheet (structural damping). See section 2.3.4 for more 
discussion about structural damping. 
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TYPE 2 PHYSICAL CONNECTING ELEMENT 
(SPACE LINK-DAMPER ELEMENT) 


V 

$ , 

$LIST3 

ITYPE=2 , 

element number 

I LEM= , 

I-end J-end 

point number point number 
JT= , , 

Area (in 2 ) 

TAREA= , 

Young's Modulus (psi) 

TYOUNG 3 , 

1 3 damping 
0=no damping 

I DAMP 3 , (0 assumed) 


If I DAMP 3 1 , complete the input for one of the following two options: 

Option 1: For damping based on Q-factor and selected frequency, enter the 

following C If QELEM is entered but not QFREQ, the frequency defined 
by the ISF variable on type A input sheet is used). 

Q-factor 
QELEM= , 

frequency (hertz) 

QFREQ= , 

Option 2: For damping based on the translational damping coefficient c along 

the axis of the link element, enter the following: 

C (lb-sec/in) 

TDRATE= , 


Z (vertical 
i direction) 


Space Link (Truss ) 


Y (horizontal 
di rection ) 


X (axial direction, positive 
forward) 
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2. 3. 4. 3 Rub Element (Type F) 


The rub element allows the mathematical modeling of the nonlinear tip rub, which 
includes dead band displacement (input variable DBAND) and a cubic nonlinear 
factor (input variable AK). The cubic nonlinear factor is a new feature that 
has been added to TETRA 2. This element can be used to model rotor to case rubs 
as well as rotor to rotor rubs. In the former case, the I end must be on the 
rotor and the J end on the case. In the latter case, the I end must be on the 
inner rotor and the J end on the outer rotor. 
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TYPE 3 PHYSICAL CONNECTING ELEMENT (RUB ELEMENT) 


Keff 


e 0 (dead band) 


~T1 

I 

Kpff - SKU-^vt) [1+AK( |a|-£ 0 )2] where 1 a | - relative deflection between 
i A l line of centers 

$ 

$LIST3 

ITYPE-3, 

I LEM* , element number 

I end J end 

point number point number 
JT* , , 


Note: If rotor-case rub, the I end must be on the rotor and the J end on the 

case. If rotor-rotor rub, the I end must be on the inner rotor and the J 
end on the outer rotor. 

DBAND* , radial dead band e 0 (mils) between rotor and case (if 

rotor-case rub) or between inner rotor and outer rotor (if 
rotor-rotor rub) (0 assumed) 

SK= , linear factor (lb/in) for radial spring constant K e ff (becomes 

active on closure) (0 assumed) 

AK= , nonlinear factor (in' 2 ) for radial spring constant K e ff (becomes 

active on closure) (0 assumed) 

CC= , damping coefficient C (lb-sec/in) (becomes active on closure) 

(0 assumed) 
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2. 3. 4. 4 Engine Support-Links Element (Type G-l through G-31 

The engine support element is a multipoint, multidirection variable geometry 
element that provides the capability to model the complex load paths between the 
engine case and the pylon or ground and thus allows for the direct modeling of 
actual engine mount structures. In addition, this element couples the flexible 
and rigid body centerline modal subsystem that represents the engine case to the 
support links that comprise the mounting system through the case flexibility. 

The case flexibility is described by the case distortion flexibility rates SKV, 
SKH, and SKA for the vertical, horizontal, and axial directions respectively. 

The flexibility rates are the reciprocals of the case spring rates that 
represent local distortion of the engine case under mount reaction loads. 
Multipliers, defined by the input variables AM, BM, CM, DM, EM, and FM, are used 
to proportion these flexibility rates as per the diagram shown on the type G-2 
input sheet. If it is desired to cut a load path at a case attach point J or K 
in a given direction, then the applicable multiplier can be set equal to zero. 
The load paths between the engine case and the pylon or ground can be defined 
with up to 8 link elements. In defining these load paths, no more than 4 
distinct points can be specified on the pylon or ground side. This means that 
more than a single link may be used to connect a case point to a pylon or ground 
point. The coordinates for the points that correspond to the 0 and K points on 
the case are defined on the Type B input sheet if these points are connected to 
links. 

The damping, if present, is defined by the input q-factor and input selected 
frequency (non-structural damping), or is defined by the input q-factor 
and either the independent rotor speed or the steady state forcing frequency as 
determined by variable ISF on the type A input sheet (structural damping). For 
more discussion on structural damping, see section 2.3.4. For more discussion 
on the type 4 element, see reference 1. 
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TYPE 4 PHYSICAL CONNECTING ELEMENT (ENGINE SUPPORT-LINKS ELEMENT) 

Some example configurations that can be modeled with the type 4 element: 



$ 

SLIST3 

ITYPE=4, 

ILEM= » element number 


I-end 

point number 
(engine centerline) 


J-end 

point number 
(Case point) 


K-end 

point number 
(Case point) 


JT= , » » 

Point I is attached to subsystem. Points J and K can connect to link load paths or 
subsystems or ground. 


$KV= , Vertical direction case distortion 

flexibility rate (in/lb) 


SKH= , Horizontal direction case distortion 

flexibility rate (in/lb) 


SKA= , Axial direction case distortion 

flexibility rate (in/lb) 
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TYPE 4 PHYSICAL CONNECTING ELEMENTS (ENGINE SUPPORT-L I NKS ELEMENT) (Continued) 


Enter multipliers a, b, c, d, e, and f to proportion the spring rates: 


V 

AM= 

BM= 

CM= 

DM= 

EM= 

FM= 


Rate 

PQ] 

NT 

J 

K 

Kv 

a 

c 

kh 

b 

d 

ka 

e 

f 


a + c = 1 .0 \ 
b + d = 1.0 < 
e + f = 1.0 J 


This restraint must be followed 
(zero values are permissible). 


Input one line for each link load path (maximum of 8). If no link load paths, 
omit this input. 


Case side point 

Other side point 



number (must 

number (no more 


Young's 

correspond to case 

than 4 distinct 


Modulus 

point J or K) 

point numbers) 

Area (in^) 

psi 


TLP(1 ,1 )= 
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TYPE 4 PHYSICAL CONNECTING ELEMENT (ENGINE SUPPORT-LINKS ELEMENT) (Continued) 

Enter the following only if you want damping based on Q-factor and selected 
frequency. If no damping desired, omit this input. If QELEM is entered but not 
QFREQ, the frequency defined by the ISF variable on type A input sheet is used. 
Q-factor 

QEL EM= , 

frequency (hertz) 

QFREQ= . 
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2. 3. 4. 5 Uncoupled Point Spring-Damper Element (Type H-l and H-21 

This input allows the connection of two points with a set of uncoupled springs 
and dampers. This element is typically used to connect the centerlines of 
beam-like modal subsystems and provides load paths in three translational and 
two rotational directions. Because of the lack of load path coupling, good 
modeling practice infers that the points being connected by this element should 
be coincident in space. 

Damping, if present, can be specified by the coefficient input (non-structural 
damping), by the input q-factor and selected frequency (non-structural damping), 
or by the input q-factor and either the independent rotor speed or the steady 
state forcing frequency as determined by variable ISF on the type A input sheet 
(structural damping). For more discussion about structural damping, see section 
2.3.4. 
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I and J points are 
coincident . 


^(vertical direction) 



(horizontal 

direction) 


(axial direc- 
tion, positive 
forward) 


I Z (vertical direction) 

^ y (horizontal direction) 

e 

x y 

(axial direction, 
positive forward) 




K (lb /in) 
x 

K y (lb/in) 

K (lb /in) 

2 ( in-lb } 
rad ; 

K 0z ( rad ; 


V 

$ 

$LIST3 
ITYPE-5 , 


ILEM= 

, element nunber 

I-end 

point 

JT* 

J-end 

number point number 

9 * 

Enter the following spring constants 

XS * 

, (0 assumed) 

YS - 

, (0 assumed) 

ZS = 

. (0 assumed) 

TYS * 

, (0 assumed) 

TZS - 

, (0 assumed) 
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TETRA 

TYPE 5 PHYSICAL CONNECTING ELEMENT 
UNCOUPLED POINT SPRING - DAMPER ELEMENT (CONTINUED) 


V 

I DAMP* 


l*damping 
0*no damping 


, (0 assumed) 


If I DAMP* 1 , complete the input for one of the following two 
options : 


Option 1: For damping based on Q-factor and selected 

frequency, enter the following (if QELEM is entered but not 
QFREQ, the frequency defined by the ISF variable on type A 
Q-factor input sheet is used). 

QELEM* , 

frequency (hertz) 

QFREQ* , 



Option 2: For damping based on damping coefficient definition, 

enter the following: 


XD = 
YD = 
ZD = 


TYD* 


TZD* 


, (0 assumed) 
, (0 assumed) 
, (0 assumed) 
, (0 assumed) 
, (0 assumed) 


4 
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2. 3. 4. 6 Squeeze Film Damper Element (Type H-3) 

This input provides the capability to model a squeeze film damper. Unlike the 
other elements, which can be used for both transient and steady state analyses, 
the squeeze film damper element can only be used for transient analysis runs. 

If this input is present for a steady state analysis run, it is simply ignored. 
For further details about the squeeze film damper element, see reference 3. 
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TYPE 6 PHYSICAL CONNECTING ELEMENT 
(SQUEEZE FILM DAMPER ELEMENT! 

Note: This input applies for transient analysis runs only. For steady state 

analysis runs, the below input is ignored. 

V 

$ 

SLIST3 

ITYPE-6, 

ILEM= , element number 

I end J end 

point number point number 
JT= , 

Note: The I end should correspond to the rotor or inner rotor. The J end 

should correspond to the stator or outer rotor. 

BD= , nominal damper annulus diameter (in) 

BL= , nominal damper annulus length (in) 

BC= , damper annulus radial clearance (in) 

VISO , damper lubricant viscosity (reyns) 

PVAP= , film rupture (vaporization) pressure (psi) 

PB1= , specified boundary pressure at port 1 (psi) 

PB2= , specified boundary pressure at port 2 (psi) 

TH1= , angular position of port 1 (degrees) 

TH2= , angular position of port 2 (degrees) 

NGRID= , number of finite difference grid points per damper arc 

(should be odd) 

NP0RT= , number of lubricant ports (1 or 2) 

NS0LN= , 1 - long bearing solution used 

2 = short bearing solution used (parabolic) 
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2.3.5. Restart. Time Integration, and Print/Plot Multiples (Type 1-1) 

First, the two lines $ and SLIST4 must be included regardless of whether the run 
is for a transient or a steady state analysis. The $ is needed to signal to the 
program the end of the previous namelist section, and $LIST4 signals the program 
the beginning of the LIST4 namelist section. 

The next few input variables on the type 1-1 input sheet apply for transient 
analysis runs only, and are ignored if the run is a steady state analysis run. 

If the run is a restart run, the user must enter ISTART*1 and, if desired, enter 
a restart time via the RTIME input variable. The restart time entered must 
correspond to one of the times for which output was printed on the initial run. 
If the run is a restart run and this variable is omitted, the program assumes 
that the restart time is the final time for which output was printed on the 
initial run (provided that the initial run did not terminate prematurely). 
Restart capability applies for transient analysis runs only. For more details 
about restart, see sections 2.1 and 3.1. 

Next, for transient analysis runs the user must input the time step (DELTA input 
variable). The time step should be made equal to about 1/40 of the smallest 
period of oscillation. Making the time step much larger than this runs the 
risk of numerical problems. 

Transient analysis runs begin at time equal 0 (unless the run is a restart run, 
in which case the run begins at the restart time), then time accumulates. The 
final time (TFINAL input variable) must be inputted so that the program knows 
when to stop. It is recommended that the user choose TFINAL such that the 
program will do a small number of time steps until the user is sure that his or 
her input is correct, so as to avoid possible costly garbage runs. 

The next input variable is the transient analysis print multiple IPRMUL. This 
variable governs the number of time steps for which output is printed and for 
which data gets written onto the output restart file. If IPRMUL * 100, then one 
out of every 100 time steps computed is printed out and written onto the restart 
file. Similarly, the transient analysis plot multiple IPLMUL governs the number 
of time steps that get written onto the output plot file. If IPLMUL = 10, then 
one out of every 10 time steps is written onto the plot file. 

Example: If DELTA*. 0001 seconds, TFINAL*. 5 seconds, IPRMUL*1000, and IPLMUL=20, 

then computations are made for 5001 time steps (starting with time * 0 and 
ending with time * .5). Output is printed for 6 time steps (0, .1, .2, .3, .4, 
and .5). Output is written to the plot file for 251 time steps (again starting 
at time = 0 and ending with time * .5). 

It is recommended that the user pick IPRMUL such that no more than 10 time steps 
are printed, in order to avoid being buried in printed output. Actually, 
printout of about 3 time steps is normally sufficient. The printed output is 
useful to verify that there are no problems, but it is the plots that are most 
important. It is a good idea to check a couple of points on each plot produced 
versus the printed output to be sure the plots are correct. 
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The user should pick IPLMUL so that there are enough points (time steps) to 
define a smooth curve when the data is plotted. Usually about 200 or 300 points 
is enough to define a smooth curve. The user might need more points than this 
for an orbit type plot, especially if TFINAL is large. However, only very 
rarely would more than 1000 points be needed. The user is cautioned against 
requesting too many points (time steps) to be written to the plot file. The 
size of the output plot file is approximately proportional to the number of time 
steps included, and requesting too many points can make for very large plot 
files. Also, the more points on a plot file, the more time it takes for the 
plotting program to do its stuff. 

The final variable on the type 1-1 input sheet is the steady state analysis 
print multiple, JPRMUL. This variable only applies for steady state analysis 
runs, and is ignored for transient analysis runs. This variable governs the 
number of forcing frequencies for which output is printed for a steady state 
analysis run, and is similar to the transient analysis print multiple IPRMUL. 

If JPRMUL = 50, then output is printed for one out of every 50 of the forcing 
frequencies for which a solution is computed. 

Note that for steady state analysis runs, all forcing frequencies are written to 
the plot file, so there is no steady state plot multiple analogous to the 
transient analysis plot multiple IPLMUL. This is because, for a transient 
analysis, one must find a solution for more time steps than are needed to 
produce a smooth curve in order to avoid numerical difficulties. However, for 
steady state runs there is no such restriction. For a steady state analysis, 
there is no reason for solving for any more forcing frequencies than are needed 
so that a smooth plot can later be obtained. 
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Note: The following two lines must be included regardless of whether the run is 

for a transient or a steady state analysis. 

W 

$ 

$LIST4 


RESTART. TIME INTEGRATION. AND PRINT/PLOT MULTIPLES 

Note: The below input applies for transient analysis runs only. For steady 

state analysis runs, the below input is ignored (a steady state run 
cannot be restarted) . 

0»new run 
l»restart run 

ISTART= , restart indicator (integer) (0 assumed) 

If the run is a restart run (ISTART=1), enter (if desired) the restart time. 
This time must correspond to one of the times for which output was printed on 
the initial run. If the run is a restart run and this input variable is 
omitted, the program assumes that the restart time is the final time for which 
output was printed on the initial run (provided that the initial run did not 
terminate prematurely). 

RTIME= , restart time (seconds) for transient analysis runs 

DELTA= , time step (seconds) for transient analysis runs 

TFINAL= , final time (seconds) for transient analysis runs 

IPRMUL= , print multiple for transient analysis runs (integer) 

IPLMUL= , plot multiple for transient analysis runs (integer) 

Note: The below input applies for steady state analysis runs only. For 

transient analysis runs, the below input is ignored. (For steady 
state analysis runs, all forcing frequencies are written to the plot 
file, so there is no steady state plot multiple). 

JPRMUL= , print multiple for steady state analysis runs (integer) 


$ (Include if this is the last line of the input file, otherwise omit). 
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2.3.6 Steady State Nonlinear Analysis Input (Type 1-2) 

Input sheet 1-2 applies only for a steady state analysis run in which at least 
one nonlinear type 3 physical connecting element (rub element) is present. For 
other runs, this input is ignored. Also, for each forcing frequency, the 
program makes a quick check to see if a rub is present (that is, the rub element 
dead band has been exceeded) for at least one of the rub elements. If no rub is 
present for any of the rub elements, then iteration is not needed to find the 
solution, and so this input is not used for that forcing frequency. See section 
4.5.4 of volume 1 for further details. 

First, the user may specify which subroutine from the IMSL library is to be used 
to solve the nonlinear equations (by iteration). If the user enters IROUT = 1 
or omits this variable, IMSL subroutine ZSCNT is used. If the user specifies 
IROUT = 2 IMSL subroutine ZSPOW is used. ZSCNT is based on the secant method 
for solving nonlinear simultaneous equations. ZSPOW is based on Powell’s hybrid 
algorithm, which is a variation of Newton’s method. The authors’ experience has 
been that both subroutines usually provide good results. Only for some cases 
with very large deflections (one inch or higher) were convergence problems 
encountered. For the few cases where convergence problems were encountered, 
subroutine ZSCNT usually worked better than ZSPOW. 

Next, input variables NSIG and ITMAX control how many iterations are performed. 
Iteration ceases if two successive approximations to a given root agree in the 
first NSIG digits. Assumed values are 3 for NSIG and 200 for ITMAX. Neither 
input variable has much influence on program efficiency. 

Finally, input variable GUESS can be used to input initial guesses for the rub 
element relative displacements (the unknowns in the iterating equations). If 
initial guesses are entered at all, they must be entered for all rub elements, 
and must be entered in the order that the rub elements appear in the input file. 
Of course, the initial guess input will only be used if the quick check that the 
program makes determines that iteration is needed, as outlined in section 4.5.4 
of volume 1. The initial guess input also only applies for the very first 
solution (forced frequency). For a detailed discussion of how initial guesses 
are arrived at for all forced frequencies, see section 4.5.5 of volume 1. 

If not entered via input variable GUESS, the initial guesses for the rub element 
relative displacements (if they are needed) will be those that would result 
assuming that there are no rub element rub element forces for any of the rub 
elements. This is usually the best place to start anyway, so the GUESS input 
will usually not be needed. Sometimes, however, the user might decide to use 
the GUESS input to enter starting values that he or she has found out about from 
previous runs. More than one solution may be possible for a given forcing 
frequency, and with knowledge gained from previous runs, the user might use the 
GUESS input to steer the program to converge to the solution desired. 
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STEADY STATE NONLINEAR ANALYSIS INPUT 

Note: The following input applies only for a steady state analysis in which at 

least one nonlinear type 3 physical connecting element (rub element) is 
present. Otherwise, the below input is ignored. 

If desired, specify the subroutine to be used to solve the nonlinear equations. 

IROUT= , 1 - IMSL subroutine ZSCNT (assumed value) 

2 - IMSL subroutine ZSPOW 

Enter the iteration stopping criteria below. Iteration stops (a root is 
accepted) if either of the following criteria is satisfied. 

NSIG= , A root is accepted if two successive approximations to a given root 

agree in the first NSIG digits (3 assumed) 

ITMAX= , Maximum number of iterations (200 assumed) 

Enter, if desired, initial guesses for the rub element relative displacements 
(the unknowns in the iterating equations). If initial guesses are entered at 
all, they must be entered for all rub elements. Enter one line (four values) 
for each rub element, in the same order that the rub elements appear in the 
input file. If not entered, the initial guesses will be those that would result 
assuming that there were no rub element forces for any of the rub elements. 
These initial guesses apply to the first solution (forced frequency) only. 


RELATIVE DISPLACEMENT COMPONENT (INCHES) 


VERTICAL 

HORIZONTAL 

VERTICAL 

HORIZONTAL 

COS 

COS 

SIN 

SIN 


GUESS= 


$ (Include if this is the last line of the input file, otherwise omit). 
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2.3.7 Rotor Speed and Rate Incut (Type J-l and J-21 

Input sheets J-l and J-2 are used for the rotor speed and rate input. This 
input is required if unbalance forces, gyroscopic forces, or nonstructural 
damping which is tied to the independent rotor speed is present. If this is not 

the case, the user can skip input sheets J-l and J-2 altogether. 

If rotor speed input is required, the user must specify which rotor is the 
independent rotor by setting IROTI to 1 or 2. Rotor 1 corresponds to subsystems 
1, 2, and 3 and rotor 2 corresponds to subsystems 4, 5, and 6. The 
independent rotor is the rotor for which the BEGRPM, BEGTIM, and TRHIS input is 
entered for a transient analysis run, or for which the BRPM, FRPM, and SRPM 
input is entered for a steady state analysis run. 

If rotor speed input is required and the run is a transient analysis run, the 

beginning time BEGTIM, beginning rpm BEGRPM, and time/rate table TRHIS must then 
be entered, which defines the speed and rate history for the independent rotor. 
There must be at least one and no more than 10 speed segments in the time/rate 
history table. However, if rotor speed input is required and the run is a 
steady state analysis run, the beginning rpm BRPM, final rpm FRPM, and step rpm 
SRPM must then be entered. 

If rotor speed input is required and another rotor is also present, the other 
rotor is referred to as the dependent rotor, because its speed is a function of 
the independent rotor speed. For the dependent rotor, the user should enter the 
coefficients A, B, C, and D relating the dependent rotor speed Y (rpm) to that 
of the independent rotor speed X (rpm), where Y * AX 3 + BX 2 + CX + D. 


59 



TETRA 


Page of 

NAMELIST 
Type J-l 


ROTOR SPEED AND RATE INPUT 

Note: This input is required if unbalance forces, gyroscopic forces, or 

nonstructural damping that is tied to the independent rotor speed is 
present. 

Enter the independent rotor number (rotor for which the following input defining 
the speed applies). Permissible values are 0 (static case -unbalance forces, 
gyroscopic forces, and nonstructural damping which is tied to the independent 
rotor speed not considered), 1 (rotor corresponding to subsystems 1, 2, and 3), 
and 2 (rotor corresponding to subsystems 4, 5, and 6). Must be integer. 

7*7 

IROT 1= , (0 assumed) 

Note: The remainder of this page applies for transient analysis runs only. For 

steady state analyses, the below input is ignored. 

If IROTI = 1 or IROTI - 2, enter the following: 

BEGTIM= , beginning time for the first speed segment (applies to the 

independent rotor) 

BEGRPM= , beginning speed (rpm) for the first speed segment (applies to 

the independent rotor) 

If IROTI * 1 or IROTI ■ 2 enter the following table in chronological order 
(applies to the independent rotor ) (maximum of 10 segments): 


Ending Time 

Rate 

(seconds) 

(rpm/sec) 


TRHIS* 


$ (Include if this is the last line of the input file, otherwise omit). 


Segment 1 
Segment 2 
Segment 3 
Segment 4 
Segment 5 
Segment 6 
Segment 7 
Segment 8 
Segment 9 
Segment 10 
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TETRA 


Page of 

NAMELIST 
Type J-2 


ROTOR SPEED AND RATE INPUT (Continued! 

Note: The following input applies for steady state analysis runs only. For 

transient analysis runs, the following input is ignored. Because (for 
steady state analysis runs only) all unbalance points must lie on the 
independent rotor, the exciting frequency will correspond to the 
independent rotor speed if unbalance loads are present. 

If IROTI * 1 or IROTI - 2, enter the following: 

V 

BRPM= , beginning rpm for the independent rotor 

FRPM= , final rpm for the independent rotor 

SRPM= , step rpm 


Note: The following input applies both for transient analysis and steady state 

analysis runs. 

If IROTI = 1 or IROTI » 2 and a second rotor is present, input the following 
coefficients relating the second (dependent) rotor speed Y (rpm) to the 
independent rotor speed X (rpm), where Y=AX 3 + BX 2 + CX + D: 

A= , (0 assumed) 

B* , (0 assumed) 

C* , (0 assumed) 

D= , (0 assumed) 


$ (Include if this is the last line of the input file, otherwise omit). 
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2.3.8 Applied Loads (Type K-l through M-2) 

2. 3. 8.1 Unbalance Loads (Type K-l and K-2) 


Input sheets K-l and K-2 are for unbalance load input. To have unbalance loads, 
rotor speed input (type J-l and J-2) is needed - otherwise, the unbalance load 
input is ignored. 

For transient analysis runs, the unbalance loads are entered via the UNBAL input 
variable. For each transient analysis unbalance load, four values must be 
entered - the time of birth (when the unbalance load becomes active), the global 
point number on the unbalanced rotor (the point must lie on rotor 1 or rotor 2), 
the magnitude of the unbalance (gm-in), and the phase angle (degrees) (see the 
figure on the type K-l input sheet). For steady state analysis runs, on the 
other hand, the unbalance loads are entered via the UBALSS input variable. Only 
three values must be entered for each steady state unbalance load - the same 
values as for a transient analysis except that the time of birth is omitted. 
Whether the run is transient or steady state, there can be from 0 to a maximum 
of 20 unbalance loads. 

The unbalance load input is quite flexible. For transient analyses, the time of 
birth can be the same for different birth events (time of birth does not apply 
for steady state analyses). Also, for both transient and steady state runs the 
same point can be referenced more than once if desired. If more than one 
unbalance load is specified for the same point, the unbalance loads are added 
together to get the total unbalance load. For example, the user could make a 
transient analysis run where a nominal unbalance load is applied at the engine 
fan location, then at a later time, after transient effects have died out, 
introduce a much larger unbalance due to a fan blade loss. 

For a transient analysis restart run (steady state runs do not have restart 
capability), any unbalance loads that were present in the original run should be 
left in for the restart run. (Otherwise, the unbalance loads would disappear 
for the new run). Additional birth events that become active after the restart 
time may be added to the UNBAL input for a restart run, however. 
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TETRA 


Pace of 

NAMELIST 
Type K-l 


* 



NOTE: To have unbalance loads, rotor speed input (type J-l and J-2) is 

required - otherwise, the unbalance load input is ignored. 

NOTE: The following input applies for transient analysis runs only. For 

steady state analysis runs, the following input is ignored. 

If unbalance loads are desired, enter the following (maximum of 20 lines): 


Time of 

Point nunber on 

Magnitude 

Phase Angle 4> 

birth (seconds) 

unbalanced rotor 

(gm-in) 

(degrees) 



If restart run and the time of birth is less than the restart time, then the 
unbalance load continues active for the restart run. Time of birth can be 
the same for different birth events. Also, the same point can be referenced 
more than once if desired. 


S (Include if this is the last line of the input file, otherwise omit). 
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TETRA 


Page of 

NAMELIST 
Type K-2 


APPLIED LOADS 

UNBALANCE LOADS (Continued) 

Note: The following applies for steady state analysis runs only. For transient 

analysis runs, the following input is ignored. Note that the variable 
name for steady state unbalance loads (UBALSS) is different than that for 
transient unbalance loads (UNBAL). Also note that all unbalance points 
must be located on the independent rotor (the independent rotor and its 
speed are specified on input sheets J-l and J-2). 

If unbalance loads are desired, enter the following (maximum of 20 lines): 


Global Point 



Number on 



Unbalanced 


Phase 

(Independent) 

Magnitude 

Angle <P 

Rotor 

gm-in 

degrees 




UBALSS* 



$ 


(Include is this is the last line of the input file, otherwise omit). 



2. 3. 8. 2 P*cos ( ait 1 and P*sin(a)tl Loads (Type L-ll 


The P*cos(o)t) and P*sin(wt) load input applies for transient analysis runs only. 
For steady state analysis runs, this input is ignored. If this type of input is 
not desired, this input sheet can be skipped altogether. If these loads are 
desired, the user must enter 5 values for each load via the CS input variable as 
shown on input sheet L-l. As for the unbalance load input, the same point can 
be referenced more than once, in which case the loads are added together to get 
the total load. There can be from 0 to a maximum of 30 loads. 

Definitions are as follows: 

P « force amplitude, (lb or in-lb) 

a) - frequency (hertz) 

t - time (seconds) 
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TETRA 


Page of 

NAME! TIT 
Type L - 1 


APPLIED LOADS 

P*CC5 (ait; and P*sin (tut) L0 £ Dr< 



Y (horizontal 
direction) 


GLOBAL 

DIRECTION 

NUMBER 

DIRECTION 

1 

Z 

2 

6y 

3 

Y 

4 

®Z 

5 

x 

6 

j ex 

i 


X (axial direction, positive forward) 


NOTE: The following applies for transient analysis runs only. For 

steady state analysis runs, the following input is ignored. 


If P*cos (wt) or P*sin (wt) loads are desired, enter the following 
(Maximum of 30 lines): 


GLOBAL 

TYPE 

AMPLITUDE 

FREQUENCY 

6L0BAL 

POINT 

(1«C0S 

P 

03 

DIRECTION 

NUMBER 

2-SIN) 

(lb or in-lb) 

(hertz) 

NUMBER 


cs= 



S (Include if this is the last line of the input file, otherwise omit). 
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2. 3. 8. 3 P*cos(^t+ji) Loads (Type L-21 


P*cos(ojt+ <f>) loads apply for steady state analysis runs only. For transient 
analysis runs, this input is ignored. If this type of loading is not desired, 
the type L-2 input sheet can be skipped altogether. If this type of loading is 
desired, 4 values are entered for each load via the COSSS input variable, as 
shown on input sheet L-2. The same point may be referenced more than once, in 
which case the loads are added together to get the total load. There can be 
from 0 to a maximum of 30 loads. 

Definitions are as follows: 

P * force amplitude (lb or in-lb) 

<*> - forcing frequency (hertz) 

t * time (seconds) 

<f> * phase angle (degrees) 
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ieim 


Page of 

NAMELIST 
Type L-2 


APPLIED LOADS 
P*cos ( tiit+d) I LOADS 

Note: The following input applies for steady state analysis runs only. For 

transient analysis runs, the below input is ignored. 

If P*cos (iot+<j>) loads are desired, input the following (maximum of 30 lines): 


Global 

Ampl itude 

Phase 

Global 

Point 

P 

Angle <J> 

Direction 

Number 

lb or in-lb 

degrees 

Number 


V 


cosss* 


If steady state analysis run and unbalance loads (input sheet K-2) are present 
in addition to P*cos(u)t+<j>) loads, the program will set the forcing frequency 
for the P*cos( (J jt+tj)) loads equal to the speed of the unbalance (independent) 
rotor, and the. below input, if present, is overridden. If, however, P*cos(-t+s) 
loads are present but not unbalance loads, you must enter the below input to 
define the forcing frequency w for the P*cos(tot+<|>) loads. A steady state 
analysis will then be performed at each of the below frequencies for every rotor 
speed defined on the type J-2 input sheet. 

BW= , beginning forcing frequency (rpm) 

FW= , final forcing frequency (rpm) 

SW= , step forcing frequency (rpm) 


$ (Include if this is the last line of the input file, otherwise omit). 
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2. 3. 8. 4 Time-Force History Loads (Type M-l and M-21 


Input sheet M-l and M-2 apply for transient analysis runs only. For steady 
state analysis runs, this input is ignored. If time-force history loads are not 
desired, input sheets M-l and M-2 can be skipped altogether. 

If time-force history loads are desired, the user must specify the global point 
* number, global direction number, and table number for each time-force history 

load via the NTF input variable on input sheet M-l (values must be integers). 

The table number must be a value between 1 and 10. There can be from 0 to a 
. maximum of 30 time-force history loads entered in the NTF array. 

Then, for each table number referenced in the NTF array, the array table is 
entered via TABLE (1,1, NT), where NT is the table number referenced in the NTF 
array. Entries in each table must be in order of increasing time. If the first 
time-force entry in the table is for a time other than zero, then the force is 
assumed to be zero up to the time of the first table entry. If the last 
time-force pair in the table is for a time less than the final time of the run, 
then for times greater than that of the last table entry the force is assumed 
constant and equal to the force for the last table entry. A table may have only 
one time-force pair entry, in which case the force is zero till the entered time 
and afterward equal to the entered force. There can be a maximum of 10 tables, 
and a maximum of 10 time-force pairs in each table. 
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TETRA 


Page of __ 

NAMELIST 
Type M-l 


APPLIED LOADS 
TIME-FORCE HISTORY LOADS 


Z (vertical direction) 



X (axial direction, positive forward) 


GLOBAL 

DIRECTION 

NUMBER 

1 

DIRECTION 

1 

Z 

2 

0 Y 

3 

Y 

4 

e z 

5 

X 

6 

0X 


Note: Sheets M-l and M-2 are for transient analysis runs only. For steady 

state analysis runs, this input is ignored. 


If time-force history loads are desired, enter the following (values must be 
integers) (maximum of 30 lines): 


GLOBAL 

GLOBAL 

TABLE NUMBER 

POINT 

DIRECTION 

NT (VALUE 

NUMBER 

NUMBER 

BETWEEN 1 & 10) 


V 

NTF = 
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TETRA 


Page of 

NAMELIST 
Type M-2 


APPLIED LOADS 

TIME-FORCE HISTORY LOADS (Continued! 

For each time-force table, enter the following. Entries must be in order of 
increasing time. Substitute the actual table number in place of the dimension 
NT below. A maximum of 10 tables and 10 time-force pairs in each table is 
permitted. 


TIME 

FORCE 

seconds 

lb or in-lb 




TABLE(1,1,NT)= 


$ 


(Include if this is the last line of the input file, otherwise omit). 




2.3.9 Gyroscopic Loads (Type N) 

Input sheet N is for the gyroscopic load input. This input applies for both 
transient and steady state analysis runs. For each gyroscopic load, the point 
number on the rotor and the polar moment of inertia (lb-in*) is entered. There 
can be a maximum of 30 gyroscopic loads. Since this input models the cross-axis 
coupling forces associated with Coriolis acceleration, both the vertical and 
horizontal subsystems must be included for the rotor(s). 


TETRA 


Page of 

NAMELIST 
Type N 


GYROSCOPIC LOADS 


Note: To have gyroscopic loading, rotor speed input (type J-l and J-2 input 

sheets) must be included. Otherwise, the gyroscopic load input is 
ignored. 

Note: For gyroscopic loading, both the vertical and horizontal subsystems must 

be included for the rotor(s). 

If gyroscopic loads are desired, enter the following (maximum of 30 lines): 


POINT 
NUMBER 
ON ROTOR 


POLAR MOMENT 
OF^INERJIA I p 



GYR0= 


$ (Include is this is the last line of the input file, otherwise omit). 
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2.3.10 Data for Plot file (Type P-1 and P-2) 

Input sheets P-1 and P-2 are used to control what data (if any) gets written to 
the output plot file. These input sheets apply for both transient and steady 
state analysis runs. 

First comes input variable IPLOT. If no output plot file is wanted, the user 
should set variable IPLOT to 0 and dispense with the rest of input sheets P-1 
and P-2. If, however, the user sets IPLOT to 1 or omits this input variable, 
then an output plot file will be produced, and the user should complete the rest 
of input sheets P-1 and P-2. 

The plot file size in many cases would be huge if all the data that conceivably 
might be desired were written onto the plot file. For this reason, the NPD 
(input sheet P-1) and NEPD (input sheet P-2) input variables were added. This 
input allows the user to specify the data that he or she wants written to the 
plot file. Be sure to include this input for all the data that you might want 
to plot. 

The NPD array allows the user to select the point and direction pairs for which 
the coordinates, displacement magnitudes, velocity magnitudes, and modal force 
magnitudes get written to the output plot file. For steady state analyses, the 
displacement phase angles, velocity phase angles, and modal force phase angles 
for these point and direction pairs are also written to the output plot file. 

If the user doesn’t want to plot values for any point and direction pairs, then 
the NPD array should be omitted. A maximum of 50 point and direction pairs may 
be specified. 

The NEPD array allows the user to select the physical connecting element number, 
global point number, and global direction number trio for which the physical 
connecting element force magnitudes are written to the output plot file. For 
steady state analyses, the physical connecting element force phase angle for 
these trios are also written to the output plot file. Again, if the user 
doesn’t want to plot values for any element, point, and direction trios, he or 
she should omit this array. A maximum of 50 element, point, and direction trios 
may be specified. 

In addition, certain other data is always automatically written to the output 
plot file (provided that an output plot file is generated). This includes the 
time (transient analysis only), forcing frequency (steady state analysis only), 
independent and dependent rotor speed, independent and dependent rotor angular 
displacement (transient analysis only), generalized displacement magnitudes 
for the flexible bladed disk modes, phase angles for the flexible bladed disk 
modes (steady state analysis only), data for the flexible bladed disk(s) (if any 
were present in the TETRA run), and data for the rub elements (if any were 
present in the TETRA run). For further details about the contents of the output 
plot file, see section 3.2. 
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TETRA 


DATA FOR PLOT FILE 


Page of 

NAMELIST 
Type P-1 


Z (vertical direction) 



Y (horizontal 
direction) 


X (axial direction, positive forward) 


GLOBAL 

DIRECTION 

NUMBER 

DIRECTION 

1 

Z 

2 

®Y 

3 

Y 

4 

®z 

5 

X 

6 

®X 


IPLOT- , 0 * no plot file 

1 ■ output plot file produced (assumed value) 


Enter the following (point, direction) pairs for which data is to be written to 
the output plot file. For both transient and steady state analysis runs this 
data includes the coordinates, displacement magnitudes, velocity magnitudes, and 
modal force magnitudes. For steady state analysis runs, this data also includes 
the displacement phase angle, velocity phase angle, and modal force phase angle. 
Values must be integers. Maximum of 50 lines. 


GLOBAL 

GLOBAL 

POINT 

DIRECTION 

NUMBER 

NUMBER 


NPD* 


$ (Includ.e is this is the last line of the input file, otherwise omit). 
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TETRA 


Page of 

NAMELIST 
Type P-2 


DATA FOR PLOT FILE (Continued) 

Enter the following (element, point, direction) trios (if any) for which data is 
to be written to the output plot file. For both transient and steady state 
analysis runs, the physical connecting element force magnitude is written to the 
output plot file for these trios. For steady state analysis runs, the physical 
connecting element force phase angle is also written to the output plot file for 
these trios. Values must be integers. Maximum of 50 lines. 


PHYSICAL 

CONNECTING 

GLOBAL 

GLOBAL 

ELEMENT 

POINT 

DIRECTION 

NUMBER 

NUMBER 

NUMBER 



NEP0= 




$ (Include if this is the last line of the input file, otherwise omit). 
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3.0 TETRA 2 OUTPUT 


3.1 Output Restart File 

If the run is for a transient analysis, an output restart file will be 
generated. This file is written to file code 24, so for transient analysis runs 
the user must have a line or two in his or her JCL (job control language) that 
assigns a name to the output restart file and associates it with file code 24. 
The output restart file thus produced may then subsequently be used as the input 
restart file for a new transient analysis restart run. 

For more information on the restart file, see section 2.1. 
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3.2 OUTPUT PLOT FILE 


Usually the user requests a plot file to be written, so that a postprocessing 
program can later be run to plot the results (plots are vital for interpreting 
the results from this program). The input variable IPLOT on input sheet P-1 
controls whether or not a plot file is produced (this default is to produce 
one) . 

The plot file is an unformatted sequential file. This plot file is written to 
file code 23. Therefore, if a plot file has been requested, the user must add a 
line or two to his or her JCL (job control language) to assign a name to the 
plot file and to associate the plot file with file code 23. 

The plot file has the potential of being a very large file, depending on how 
many time steps (for a transient analysis run) or how many forcing frequencies 
(for a steady state analysis run) are written to it. One means by which the 
size of the plot file is limited is that the user specifies (via input sheets 
P-1 and P-2) what data is to be written to this file. For more details about 
this, see the discussion of the type P-1 and P-2 input sheets. 

For a transient analysis run, the user can also limit the size of the plot file 
by controlling how many of the time steps analyzed are written to the plot file 
via the input plot multiple IPLMUL on the type 1-1 input sheet. See the 
discussion of the type 1-1 input sheet for advice on what value to use for 
IPLMUL. For a steady state analysis run, on the other hand, all of the forcing 
frequencies analyzed are written to the plot file (provided, of course, that the 
user requested a plot file). The reason that a plot multiple is provided for a 
transient analysis only is that for a transient analysis the user must run more 
time steps than needed for plotting to avoid numerical problems, whereas there 
is no such requirement for a steady state analysis. In other words, for a 
steady state analysis, there is no good reason for analyzing any more forcing 
frequencies than needed for the plot file so that a smooth plot can be later 
obtained. 

The plot file contents may be divided into sections, of which there are three 
types as follows. The first type of section is the plot file preface (see table 
3-1), which comes first and contains information telling the plot program what 
type of run (transient or steady state) and what information is present on the 
plot file (which points, directions, forces, etc.). This section is present for 
both transient and steady state analyses. This section occurs only once, and 
every line is a record for this section. 

The second section is for the transient analysis plot data (see table 3-2), and 
is present only if the TETRA run that generated the plot file was a transient 
analysis run. This section is repeated for each time step written to the plot 
file. The whole section (that is, each time step) is a record. 

Finally, the third section is for the steady state analysis plot data (see table 
3-3), and is present only if the TETRA run that generated the plot file was a 
steady state analysis run. This section is repeated for each forcing frequency 


78 



that is written to the plot file. The whole section (that is, each forcing 
frequency) is a record. 

Definitions of the plot file variables contained in all three sections (arranged 
alphabetically) are in table 3-4. 

Note that the first two types of sections (the plot file preface and the 
transient analysis plot data) are almost identical to what they were in the 
original TETRA program (which could only perform transient analyses). The only 
difference to the plot file preface is the addition of the variables NTYPE, 
IFBDR, and TINT which are tacked onto the end of the first line of the plot file 
preface. The only difference for the transient analysis plot file section is the 
addition of the variables DVM and DHM (which are useful in making orbit type 
plots). Thus, plot programs which were written for the original TETRA program 
(which was for transient analyses only) will have to be modified only slightly 
to work for transient analysis TETRA 2 runs, but will have to be modified much 
more extensively if they are also to work for steady state analysis TETRA 2 
runs. 
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Table 3-1. Plot file preface 


NP.NL1, NL2 , NRE , NEL , NPLOTP , NTYPE , IFBDR .TINT 
IP PLOT j , IDPLOTj ,XPTi , YPTj , ZPTj 


I PPL0T N p , IDPLOT md , XPTmd , YPTmd , ZPTmd 
LFBO l x , RFB01 x , AFBD1 j 

LFBD2mi 9 > RFBD2mi2 » AFBD2mi o 
I LEM3 ! , 1 LEM3 \ , I L tM3 j 

ILEM3 NR r,ILEM3 NRE ,ILEM3 NRE 
I LEM j , I PT j , IDIRj 


ilem NE l» ipt nel» idir nel 
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Table 3-2. Plot file data for a transient analysis run 


TIME , RPMI , RPMD , THETA I , THETAD 
PMAG1 , QMAG1 , PMAG2 , QMAG2 
XMAGi , VMAGj , FMODMj 


XMAG N p,VMAG N p,FMODM N p 

PS 1 1 ! , UMAG1 ! , VMAG1 j , §1 FBD1 j , S2FBD1 j , S3FBD1 j 


PSliwiij UMAG 1 
PSI2 1 ,0MAG2 1 , 




PS I 2 N i 2 > UMAG2 NL2 , VMAG2 NL2 , SI FBD2 NL2 , S2FBD2 NL2 , S3 FBD2 NL2 
RDi , CLi , Fj , DVMi , DHMi 


RDnre > CL nre , f nre , dvm nre , dhm nre 
FMAGj 


f mag N el 
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Table 3-3. Plot file data for a steady state analysis run 


WFF,RPMI,RPMD 

PMAG1 , PPH1 , QMAG1 , QPH1 

PMAG2 , PPH2 , QMAG2 , QPH2 

XMAG! , X PHASE! , VMAGi , VPHASEi , FMODMi , FMODPHi 


XMAG N p , XPHASE N p , VMAG NP , VPHASE N p , FMODM N p , FMODPH np 
UMAGP1 1 ,VMAGP1 1 ,SIG1P1 1 ,SIG2P1 1 ,SIG3P1 1 



UMAGQ1 NL1 , VMAGQ1 NL1 ,SIG1Q1 NL1 ,SIG2Q1 NL ! ,SIG3Q1 nl1 
UMAGP2 1 ,VMAGP2 1 ,SIG1P2 1 ,SIG2P2 1 ,SIG3P2 1 


UMAGP2 NL 2 ,VMAGP2 N i2,SIG1P2 NL 2,SIG2P2 NL 2,SIG3P2 NL2 
UMAGQ2 1 ,VMAGQ2 1 ,SIG1Q2 1 ,SIG2Q2 1 ,SIG3Q2 1 


UMAGQ2 NL 2,VMAGQ2 NL 2,SIG1Q2 NL2 ,SIG2Q2 NL 2,SIG3Q2 NL 2 
RX 1 ,RN 1 ,CX 1 ,CN 1 ,FX 1 ,FN 1 , DVMi , DVPi , DHMi , DHPj 


RX N re , RN N rf , cx NRE , cn nre , fx nre , fn nre , dvm nre , dvp nre , dhm nre , dhp nre 

FMAGiFPHASEi 


FMAG NEL ,FPHASE NEL 
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Table 3-4 Plot file variable definitions (alphabetical order) 


4 


VARIABLE DEFINITION 


AFBDl-j 

AFBD2-f 

CLi 

CNi 

CXi 

DHMi 

DHPj 

DVM-j 

DVP,- 

Fl 

FMAGi 

FMODMi 

FMODPH-j 

FN i 

FPHASE-j 

FXi 

IDIRi 

IDPLOT i 


Axial coordinate for the i’th local point on flexible bladed disk 
number 1. Applies for both transient and steady state analyses. 

Axial coordinate for the i’th local point on flexible bladed disk 
number 2. Applies for both transient and steady state analyses. 
Clearance (inches) for the i’th type 3 physical connecting element 
(rub element). Applies for a transient analysis run only. 

Minimum clearance (inches) for the i’th type 3 physical connecting 
element (rub element). Applies for a steady state analysis run only. 
Maximum clearance (inches) for the i’th type 3 physical connecting 
element (rub element). Applies for a steady state analysis run only. 
Magnitude (inches) of the difference between the horizontal direction 
displacement at joint I minus that of joint J for the I’th rub 
element. Applies for both transient and steady state analysis runs. 
Phase angle (degrees) for the difference between the horizontal 
direction displacement at joint I minus that of joint J for the i’th 
rub element. Applies for steady state analysis runs only. 

Magnitude (inches) of the difference between the vertical direction 
displacement at joint I minus that of joint J for the i’th rub 
element. Applies for both transient and steady state analysis runs. 
Phase angle (degrees) for the difference between the vertical 
direction displacement at joint I minus that of joint J for the i’th 
rub element. Applies for steady state analysis runs only. 

Force magnitude (pounds) for the i-th type 3 physical connecting 
element (rub element). Applies for a transient analysis run only. 
Force magnitude (pounds) for the i’th (element, point, direction) trio. 
Applies for both transient and steady state analysis runs. 

Modal force magnitude (lb or in-lb) for the i’th (point, direction) 
pair- Applies for both transient and steady state analyses. 

Modal force phase angle (degrees) for the i’th (point, direction) 
pair. Applies for steady state analysis runs only. 

Minimum harmonically averaged force magnitude (pounds) for the i’th 
type 3 physical connecting element. Applies for a steady state 
analysis run only. 

Phase angle (degrees) for the i’th (element, point, direction) trio. 
Applies for steady state analysis runs only. 

Maximum harmonically averaged force magnitude (pounds) for the i’th 
type 3 physical connecting element. Applies for a steady state 
analysis run only. 

Global direction number for the i’th (element, point, direction) trio 
for which data is written to the plot file. Applies for both 
transient and steady state analysis runs. 

Global direction number for the i’th (point, direction) pair for which 
data is written to the plot file. Applies for both transient and 
steady state analyses. 
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IFBDR 

ILEM-j 

ILEM3-j 

IPPLOTi 

IPTi 

LFBDl-j 

LFBD2j 

NEL 

NL1 

NL2 

NP 

NPLOTP 


NRE 

NTYPE 

PMA61 

PMAG2 

PPH1 

PPH2 


Indicator telling which rotor contains the flexible bladed disks. 
Permissible values are 0 (no flexible bladed disks), 1 (flexible 
bladed disks are on the independent rotor), and 2 (flexible bladed 
disks are on the dependent rotor). 

Element number for the i’th (element, point, direction) trio for which 
data is written to the plot file. Applies for both transient and 
s t ea( jy snsilysGS . 

Element number for the i’th type 3 physical connecting element (rub 
element). Applies for both transient and steady state analyses. 
Global point number for the i’th (point, direction) pair for which 
data is written to the plot file. Applies for both transient and 
s teady stcitG snElysss. 

Global point number for the i’th (element, point, direction) trio for 
which data is written to the plot file. Applies for both transient 
and steady state analysis runs. 

Local point number for the i’th local point on flexible bladed disk 
number 1. Applies for both transient and steady state analyses. 

Local point number for the i’th local point on flexible bladed disk 
number 2. Applies for both transient and steady state analyses. 
Number of (element, point, direction) trios for which data is written 
to the plot file. Applicable for both transient and steady state 
analysis runs. 

Number of local points on flexible bladed disk number 1. Applicable 
for both transient and steady state analysis runs. 

Number of local points on flexible bladed disk number 2. Applicable 
for both transient and steady state analysis runs. 

Number of (point, direction) pairs for which data is written to the 
plot file. Applicable for both transient and steady state analysis 
runs. 

Number of plot points that should have been written to the plot file. 
For a transient analysis, this equals the number of time steps that 
should have been written to the plot file. For a steady state 
analysis, this equals the number of forcing frequencies that should 
have been written to the plot file. If for any reason the TETRA run 
terminated prematurely, the actual number of plot points might be 
less than this number. 

Number of type 3 physical connecting elements (rub elements). 
Applicable for both transient and steady state analysis runs. 

Analysis type code, where NTYPE-0 signifies a transient analysis run, 
and NTYPE*1 signifies a steady state analysis run. 

Generalized displacement magnitude for the horizontal nodal diameter 
mode for flexible bladed disk number 1. Applies for both transient 
and steady state analyses. 

Generalized displacement magnitude for the horizontal nodal diameter 
mode for flexible bladed disk number 2. Applies for both transient 
and steady state analyses. 

Phase angle (degrees) for the generalized displacement for the 
horizontal nodal diameter mode for flexible bladed disk number 1. 
Applies for a steady state analysis run only. 

Phase angle (degrees) for the generalized displacement for the 
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PSIli 

PSI2-} 

QMAG.l 

QMAG2 

QPH1 

QPH2 

RD-j 

RFBDl-j 

RFBD2i 

RN-j 

RPMD 

RPMI 

RX-j 

SlFBDl-j 
SlFBD2i 
S2FBD1 i 
S2FBD2 -j 
S3FBD1 i 
S3FBD2 i 
SIGlPl-j 


horizontal nodal diameter mode for flexible bladed disk number 2. 
Applies for a steady state analysis run only. 

Angle psi (degrees) for the i’th local point on flexible bladed disk 
number 1. Applies for a transient analysis run only. 

Angle psi (degrees) for the i’th local point on flexible bladed disk 
number 2. Applies for a transient analysis run only. 

Generalized displacement magnitude for the vertical nodal diameter 
mode for flexible bladed disk number 1. Applies for both transient 
and steady state analyses. 

Generalized displacement magnitude for the vertical nodal diameter 
mode for flexible bladed disk number 2. Applies for both transient 
and steady state analyses. 

Phase angle (degrees) for the generalized displacement for the 
vertical nodal diameter mode for flexible bladed disk number 1. 
Applies for a steady state analysis run only. 

Phase angle (degrees) for the generalized displacement for the 
vertical nodal diameter mode for flexible bladed disk number 2. 
Applies for a steady state analysis run only. 

Relative displacement magnitude for the i’th type 3 physical 
connecting element (rub element). Applies for a transient analysis 
run only. 

Radial coordinate for the i’th local point on flexible bladed disk 
number 1. Applies for both transient and steady state analysis runs 
Radial coordinate for the i’th local point on flexible bladed disk 
number 2. Applies for both transient and steady state analysis runs 
Minimum relative displacement magnitude for the i’th type 3 physical 
connecting element (rub element). Applies for a steady state 
analysis run only. 

Speed (rpm) of the dependent rotor. Applicable for both transient 
and steady state analyses. 

Speed (rpm) of the independent rotor. Applicable for both transient 
and steady state analyses. 

Maximum relative displacement magnitude for the i’th type 3 physical 
connecting element (rub element). Applies for a steady state 
analysis run only. 

SIGMA1 stress (psi) for the i’th local point on flexible bladed disk 
number 1. Applies for a transient analysis run only. 

SIGMA1 stress (psi) for the i’th local point on flexible bladed disk 
number 2. Applies for a transient analysis run only. 

SIGMA2 stress (psi) for the i’th local point on flexible bladed disk 
number 1. Applies for a transient analysis run only. 

SIGMA2 stress (psi) for the i’th local point on flexible bladed disk 
number 2. Applies for a transient analysis run only. 

SIGMA3 stress (psi) for the i’th local point on flexible bladed disk 
number 1. Applies for a transient analysis run only. 

SIGMA3 stress (psi) for the i’th local point on flexible bladed disk 
number 2. Applies for a transient analysis run only. 

SIGMA1 stress (psi) for the horizontal nodal diameter mode for the 
i’th local point on flexible bladed disk number 1. Applies for a 
steady state analysis run only. 
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SIG1P2-J 

SIGIQl-j 

SIGlQ2-j 

SIG2P1 

S IG2P2 -j 

SIG2Q1 i 

SIG2Q2-J 

SIG3P1-J 

SIG3P2-j 

SIG3Q1 -j 

S IG3Q2 i 

THETAD 

THETA I 

TIME 

TINT 

UMAGl-j 

UMAG2 1 

UMAGPl-j 

UMAGP2i 


SIGMA1 stress (psi) for the horizontal nodal diameter mode for the 
i’th local point on flexible bladed disk number 2. Applies for a 
steady state analysis run only. 

SIGMA1 stress (psi) for the vertical nodal diameter mode for the i’th 
local point on flexible bladed disk number 1. Applies for a steady 
state analysis run only. 

SIGMA1 stress (psi) for the vertical nodal diameter mode for the i’th 
local point on flexible bladed disk number 2. Applies for a steady 
state analysis run only. 

SIGMA2 stress (psi) for the horizontal nodal diameter mode for the 
i’th local point on flexible bladed disk number 1. Applies for a 
steady state analysis run only. 

SIGMA2 stress (psi) for the horizontal nodal diameter mode for the 
i’th local point on flexible bladed disk number 2. Applies for a 
steady state analysis run only. 

SIGMA2 stress (psi) for the vertical nodal diameter mode for the i’th 
local point on flexible bladed disk number 1. Applies for a steady 
state analysis run only. 

SIGMA2 stress (psi) for the vertical nodal diameter mode for the i’th 
local point on flexible bladed disk number 2. Applies for a steady 
state analysis run only. 

SIGMA3 stress (psi) for the horizontal nodal diameter mode for the 
i’th local point on flexible bladed disk number 1. Applies for a 
steady state analysis run only. 

SIGMA3 stress (psi) for the horizontal nodal diameter mode for the 
i’th local point on flexible bladed disk number 2. Applies for a 
steady state analysis run only. 

SIGMA3 stress (psi) for the vertical nodal diameter mode for the i’th 
local point on flexible bladed disk number 1. Applies for a steady 
state analysis run only. 

SIGMA3 stress (psi) for the vertical nodal diameter mode for the i’th 
local point on flexible bladed disk number 2. Applies for a steady 
state analysis run only. 

Angular displacement (revolutions) for the dependent rotor. 

Applicable for transient analyses only. 

Angular displacement (revolutions) for the independent rotor. 
Applicable for transient analyses only. 

Time (seconds). Applicable for a transient analysis only. 

Time interval (seconds) between the time steps on the plot file. 

This variable applies for a transient analysis run onl,y. For a 
steady state analysis, this variable is present on the plot file but 
equals 0. 

Displacement u (inches) for the i’th local point on flexible bladed 
disk number 1. Applies for transient analysis runs only. 

Displacement u (inches) for the i’th local point on flexible bladed 
disk number 2. Applies for transient analysis runs only. 

Displacement magnitude u (inches) for the horizontal nodal diameter 
mode for the i’th point on flexible bladed disk number 1. Applies 
for steady state analysis runs only. 

Displacement magnitude u (inches) for the horizontal nodal diameter 
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UMAGQ 1 i 

UMAGQ2i 

VMAG-j 

VMAGli 

VMAG2i 

VMAGPli 

VMAGP2-j 

VMAGQ1 1 * 

VMAGQ2j 

VPHASEi 

WFF 

XMAGi 

XPHASEi 

XPTi 

YPT-j 

ZPT-j 


mode for the i’th point on flexible bladed disk number 2. Applies 
for steady state analysis runs only. 

Displacement magnitude u (inches) for the vertical nodal diameter 
mode for the i’th point on flexible bladed disk number 1. Applies 
for steady state analysis runs only. 

Displacement magnitude u (inches) for the vertical nodal diameter 
mode for the i’th point on flexible bladed disk number 2. Applies 
for steady state analysis runs only. 

Velocity magnitude (in/sec or rad/sec) for the i’th (point, direction) 
pair. Applies for both transient and steady state analyses. 
Displacement v (inches) for the i’th local point on flexible bladed 
disk number 1. Applies for transient analysis runs only. 

Displacement v (inches) for the i’th local point on flexible bladed 
disk number 2. Applies for transient analysis runs only. 

Displacement magnitude v (inches) for the horizontal nodal diameter 
mode for the i’th point on flexible bladed disk number 1. Applies 
for steady state analysis runs only. 

Displacement magnitude v (inches) for the horizontal nodal diameter 
mode for the i’th point on flexible bladed disk number 2. Applies 
for steady state analysis runs only. 

Displacement magnitude v (inches) for the vertical nodal diameter 
mode for the i’th point on flexible bladed disk number 1. Applies 
for steady state analysis runs only. 

Displacement magnitude v (inches) for the vertical nodal diameter 
mode for the i’th point on flexible bladed disk number 2. Applies 
for steady state analysis runs only. 

Velocity phase angle (degrees) for the i’th (point, direction) pair. 
Applies for a steady state analysis only. 

Steady state forcing frequency (rpm). Applicable for a steady state 
analysis only. 

Displacement magnitude (inches or radians) for the i’th 

(point, direction) pair. Applies for both transient and steady state 

analyses. 

Displacement phase angle (degrees) for the i’th (point, direction) 
pair. Applies for steady state analyses only. 

X coordinate for the i’th (point, direction) pair for which data is 
written to the plot file. Applies for both transient and steady 
state analyses. 

Y coordinate for the i’th (point, direction) pair for which data is 
written to the plot file. Applies for both transient and steady 
state analyses. 

Z coordinate for the i’th (point, direction) pair for which data is 
written to the plot file. Applies for both transient and steady 
state analyses. 
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3.3 Printed Output 


For every TETRA run, printed output is produced. The file code for the printed 
output is the standard 06 file. 

The size of the printed output depends on the value inputted for the printout 
option variable IOUT on type A input sheet. See the discussion of the type A 
input for advice on what value to use for IOUT. For a transient analysis, the 
size of the printed output also depends greatly on the value inputted for the 
transient analysis print multiple IPRMUL on type 1-1 input sheet. Likewise, for 
a steady state analysis, the size of the printed output depends greatly on the 
value inputted for the steady state analysis print multiple JPRMUL, also on type 
1-1 input sheet. See the discussion of the type 1-1 input sheet for advice on 
what values to use for IPRMUL and JPRMUL. 

For a sample printed output listing, see section 5.3.1. This section contains a 
listing of the printed output for the demonstator model with the standard form 
of the printout (I0UT=1 on type A input sheet) requested. 


4.0 PROGRAM STRUCTURE AND SUBROUTINES 


TETRA 2 represents a substantial growth in program size as well as program 
capability. The previous version of TETRA (including flexible bladed disk 
capability but with only transient analysis capability and not steady state 
analysis capability) consisted of a main routine plus 36 subroutines. TETRA 2 
consists of a main routine plus 87 subroutines. Not included in this total are 
5 IMSL library subroutines which are called directly from TETRA subroutines, and 
several additional IMSL library subroutines which are not called directly from 
TETRA subroutines, but rather from the other IMSL library subroutines in TETRA. 
In terms of lines of coding, the previous version of TETRA had about 8000 lines, 
while TETRA 2 has about 30000 lines. 

The large increase in program size is largely because the steady state analysis 
portion of TETRA proved to be more complex and requiring more steps than the 
transient analysis portion of TETRA. In particular, a large amount of coding 
was needed for the assembly of the large global matrices, which are required for 
steady state analyses but not for transient analyses. The coding for the 
nonlinear rub element capability also proved to be much more extensive for the 
steady state analysis portion than for the transient analysis portion, with a 
lot of coding needed to assemble the iterating equations and perform harmonic 
averaging. Also contributing to the increase in program size was additional 
capability added to TETRA 2, such as expanding the rub element to include a 
nonlinear cubic factor, addition of structural damping capability, added 
printout options, etc. Contributing to the large increase in the number of 
lines is the fact that more extensive documentation for each subroutine is 
included in the coding (primarily at the beginning of each subroutine). 

TETRA 2 also requires appreciably more memory than the previous version of 
TETRA. A major reason that it needs more memory is for the storage of the large 
global matrices. These global matrices are needed for the steady state analysis 
portion of TETRA 2, but were not needed in the previous version of TETRA which 
was limited to transient analyses. However, the amount of memory required by 
TETRA 2 (slightly over half a million words) is rather small compared to many 
structural analysis computer codes, such as those used for finite element 
analysis. 

The structure of the TETRA 2 computer program is shown in figures 4-1 through 
4-7. As is seen from figure 4-1, the main routine of TETRA calls four different 
subroutines: PR0CES, TILOOP, FRLOOP, and LISTPF. Three of these subroutines, 
PROCESS, TILOOP, and FRLOOP, are the control routines for the three major 
branches of the program. The branch headed up by the PROCESS subroutine 
processes the input data, which is largely common to both transient and steady 
state analyses. The branch headed by the TILOOP subroutine is for the time 
integration loop, which is the solution phase for transient analysis runs. 
Finally, the branch headed up by the FRLOOP subroutine is for the frequency 
response loop, which is the solution phase for steady state analysis runs. 
Figures 4-2 through 4-7 show the branches of the program in detail. Included 
are all the TETRA 2 subroutines and all the IMSL library subroutines that are 
called directly from TETRA 2 subroutines. Omitted for clarity are those IMSL 


89 



subroutines that are not called directly from TETRA 2 subroutines but rather 
from the other IMSL library subroutines in TETRA 2. 

Table 4-1 gives a brief description of each of the subroutines of TETRA 2, 
arranged alphabetically. Included are all the TETRA 2 subroutines and all the 
IMSL library subroutines that are called directly from TETRA 2 subroutines. 
Again, those IMSL subroutines that are not called directly from TETRA 2 
subroutines, but rather from the other IMSL library subroutines in TETRA 2 are 
omitted. 
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Figure 4-1. TETRA 2 Structure Chart 
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LISTPF - Provided that a plot file was produced and that the long form of the output was 
requested, lists at least part of the output plot file. 




Figure 4-2. TETRA 2 Structure Chart (Continued 
(Parent chart is figure 4-1) 



denotes IMSL library subroutine 

























Figure 4-3. TETRA 2 Structure Chart (Continued 
(Parent chart is figure 4-2) 
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denotes IMSL library subroutines 





Figure 4-4. TETRA 2 Structure Chart (Continued) 
(Parent chart is figure 4-2) 
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Figure 4-5. TETRA 2 Structure Chart (Continued 
(Parent chart is figure 4-1) 
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Figure 4-6. TETRA 2 Subroutine Chart (Continued) 
(Parent chart is Figure 4-1) 
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denotes IMSL library subroutine 











Figure 4-7. TETRA 2 Subroutine Chart (Continued 
(Parent chart is figure 4-6) 
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denotes IMSL subroutines 












Table 4-1. TETRA 2 Subroutines 


APFOR Calculates the applied forces for transient analyses 

APF0R6 Calculates the applied generalized forces for steady state analyses 

APFORP Calculates the applied physical forces for steady state analyses 
BACKS Back substitutes the relative displacement components found by 

iteration into the iterating equations, to provide a measure of the 
quality of the solution (steady state analyses only) 

CHECK Checks if a rub is present for any of the physical connecting elements 
and acts accordingly (steady state analyses only) 

CONEL Control routine for the branch of the program which processes the input 
data for the physical connecting elements (transient and steady state 
analyses) 

CURRT Calculates the magnitudes of the physical displacements, physical 
velocities, and modal forces for transient analyses 
CURRTS Calculates the magnitudes and phase angles of the physical 

displacements, physical velocities, and modal forces for steady state 
analyses 

ELI Computes the contributions of the current type 1 physical connecting 

element (general spring-damper element) to the global stiffness matrix, 
the global velocity matrix, and the structural damping contributions 
matrix (steady state analyses only) 

EL2 Computes the contributions of the current type 2 physical connecting 

element (space link-damper element) to the global stiffness matrix, the 
global velocity matrix, and the structural damping contributions matrix 
(steady state analyses only) 

EL4 Computes the contributions of the current type 4 physical connecting 
element (engine support-links element) to the global stiffness matrix, 
the global velocity matrix, and the structural damping contributions 
matrix (steady state analyses only) 

EL5 Computes the contributions of the current type 5 physical connecting 

element (uncoupled poiat spring-damper element) to the global stiffness 
matrix, the global velocity matrix, and the structural damping 
contributions matrix (steady state analyses only) 

ELEM1 Processes the input data for the current type 1 physical connecting 
element (general spring-damper element) (transient and steady state 
analyses) 

ELEM2 Processes the input data for the current type 2 physical connecting 
element (space link-damper element) (transient and steady state 
analyses) 

ELEM3 Processes the input data for the. current type 3 physical connecting 
element (rub element) (transient and steady state analyses) 

ELEM4 Processes the input data for the current type 4 physical connecting 
element (engine support-links element) (transient and steady state 
analyses) 

ELEM5 Processes the input data for the current type 5 physical connecting 
element (uncoupled point spring-damper element) (transient and steady 
state analyses) 
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ELEM6 Processes the input data for the current type 6 physical connecting 
element (squeeze film damper element) (transient analyses only) 

ELEMV Finds the contributions of the linear (types 1, 2, 4, and 5) physical 
connecting elements (if any) and the gyroscopic elements (if any) to 
the upper right quadrant of the global solution matrix (steady state 
analyses only). 

ELG Computes the contribution of the current gyroscopic element to the 

gyroscopic contributions matrix for the independent rotor (if the gyro 
element lies on the independent rotor) or to the gyroscopic element 
contributions matrix for the dependent rotor (if the gyro element lies 
on the dependent rotor) (steady state analyses only) 

FBD Processes data for the flexible bladed disks (if any) (transient and 
steady state analyses) 

FBDDS Calculates the magnitudes of the physical displacements and stresses of 
the local points on the flexible bladed disk(s) (if any) for transient 
analyses 

FBDDSS Calculates the magnitudes of the physical displacements and stresses of 
the local points on the flexible bladed disk(s) (if any) for steady 
state analyses 

FBDL Processes the input data for the local points on the flexible bladed 
disk(s) (if any) (transient and steady state analyses) 

FBDMAS Computes the non-diagonal terms of the global mass matrix that result 
from the flexible bladed disk(s) (steady state analyses) 

FBOSS Processes the input data for the flexible bladed disk subsystems (if 
any) (transient and steady state analyses) 

FBDVEL Computes the non-diagonal terms of the global velocity matrix that 

result from the flexible bladed disk(s) (if any) (steady state analyses 
only) 

FCN Computes the values of the iterating equations for steady state 
nonlinear analyses in which a rub has occurred 
FLEX Processes the input mode shape data for flexible subsystems (transient 
and steady state analyses) 

FLEXSS Processes the input data for the flexible subsystems (transient and 
steady state analyses) 

FORCE Calculates the forces that the physical connecting elements and 

gyroscopic elements exert on the subsystems or ground for transient 
analyses 

F0RCE1 Calculates the forces that the type 1 physical connecting elements 

(general spring-damper elements) exert on the subsystems or ground for 
transient analyses 

F0RCE2 Calculates the forces that the type 2 physical connecting elements 
(space link-damper elements) exert on the subsystems or ground for 
transient analyses 

F0RCE3 Calculates the forces that the type 3 physical connecting elements (rub 
elements) exert on the subsystems for transient analyses 
F0RCE4 Calculates the forces that the type 4 physical connecting elements 

(engine support-links elements) exert on the subsystems or ground for 
transient analyses 

FORCES Calculates the forces that the type 5 physical connecting elements 
(uncoupled point spring-damper elements) exert on the subsystems or 
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ground for transient analyses 

F0RCE6 Calculates the forces that the type 6 physical connecting elements 

(squeeze film damper elements) exert on the subsystems or ground for 
transient analyses 

FORCEG Calculates the forces that the gyroscopic elements exert on the rotors 
for transient analyses 

FORHIS Processes the input data for the time-force history loads (transient 
analyses only) 

FORSS Calculates the forces that the physical connecting elements or 

gyroscopic elements exert on the modal subsystems or ground for steady 
state analyses 

F0RSS1 Calculates the forces that the type 1 physical connecting elements 

(general spring-damper elements) exert on the subsystems or ground for 
steady state analyses 

F0RSS2 Calculates the forces that the type 2 physical connecting elements 
(space link-damper elements) exert on the subsystems or ground for 
steady state analyses 

F0RSS3 Calculates the forces that the type 3 physical connecting elements (rub 
elements) exert on the subsystems for steady state analyses 

F0RSS4 Calculates the forces that the type 4 physical connecting elements 

(engine support-links elements) exert on the subsystems or ground for 
steady state analyses 

FORSS5 Calculates the forces that the type 5 physical connecting elements 
(uncoupled point spring-damper elements) exert on the subsystems or 
ground for steady state analyses 

FORSSG Calculates the forces that the gyroscopic elements exert on the 
rotors for steady state analyses 

FRLOOP Loops over the forcing frequencies and finds the solution for each 

forcing frequency in turn. This subroutine is the control routine for 
the solution phase of steady state analyses. 

GEN Calculates the generalized forces for transient analyses 

GENDIS Calculates the generalized displacements for transient analyses 

GL0B1 Prints out (provided that the user requested the long form of the 

printout) the global stiffness matrix, the global velocity matrix, the 
structural damping contributions matrix, the independent rotor 
gyroscopic contributions matrix, and the dependent rotor gyroscopic 
contributions matrix (steady state analyses only) 

GL0B2 Completes the assembly of the global solution matrix, including the 
terms that vary with rotor speed of forcing frequency (steady state 
analyses only) 

GYROE Processes the gyroscopic load input (transient and steady state 
analyses) 

INIT Performs certain initial tasks, namely: initializing, then reading the 

LIST1 namelist input, printing out the title page, listing the input 
file while at the same time counting the number of subsystems NSUB and 
the number of physical connecting elements NELEK, listing the program 
limits, processing the data for points not on the modal subsystems, and 
(for steady state analyses) initializing several matrices to 0 
(transient and steady state analyses) 

LIMITS Prints the current program. 1 i mi ts for the input variables (transient 
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and steady state analyses) 

LINSOL Solves for the generalized displacements for a linear steady state 
analysis (one for which there are no nonlinear physical connecting 
elements) 

LIST4N Processes the LIST4 namelist input (miscellaneous input that comes 
after the LIST1, LIST2, and LIST3 namelist input sections) 

LISTPF Lists at least a portion of the output plot file (provided that an 

output plot file was produced and that the user requested the long form 
of the printed output) 

MASSM Computes the mass matrix for the rotor on which the flexible bladed 
disk(s) are located and for the flexible bladed disk(s) (transient 
analyses only) 

MODES Finds the mode shape given the modal subsystem number, local mode 

number, local point number, and local direction number (transient and 
steady state analyses) 

NLDISG Calculates the generalized displacements for steady state analyses in 
which at least one nonlinear (rub) element is present 

NLFORG Calculates the generalized forces for the type 3 physical connecting 
elements (rub elements) for steady state nonlinear analyses in which a 
rub has occurred 

NLFORP Calculates the physical forces at joint I for each type 3 physical 

connecting element (rub element) for steady state analyses in which a 
rub has occurred 

NLSOL Solves for the generalized displacements for steady state analyses that 
include nonlinear type 3 physical connecting elements (rub elements) 

PARAM1 Finds the S parameters for steady state analyses in which at least one 
type 3 physical connecting element (rub element) is present. These 
parameters are needed to check if a rub has occurred, and, if a rub has 
occurred, then these parameters are also needed for use in the 
iterating equations. 

PARAM2 Finds the A, B, C, and D parameters for steady state analyses in which 

a rub has occurred. These parameters are needed for use in the 

iterating equations. 

PHYSP Processes input data for the physical points not located on the modal 
subsystems (transient and steady state analyses) 

PLOTD Processes the input plot data and, if a plot file was requested, writes 
the initial plot file data onto the output plot file (transient and 
steady state analyses) 

PROCES Processes the input data (transient and steady state analyses) 

RBODSS Processes the input data for the rigid body modal subsystems (transient 

and steady state analyses) 

RBODY Processes the mode shape input data for the rigid body modal subsystems 
(transient and steady state analyses) 

ROPROP Calculates the rotor properties (rotor speed, acceleration, and angular 
displacement for the independent and dependent rotors) (transient 
analyses only) 

SCAN Scans over the subsystems, points, directions, and modes and saves data 
needed later to calculate the generalized forces (transient and steady 
state analyses) and the physical displacements, physical velocities, 
and modal forces (transient analyses only) 
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SINCOS Processes the input data for the transient analysis P*cos(wt) and 
P*sin(wt) loads and the steady state analysis P*cos(wt+PHI) loads 

SOLVE Solves the nonlinear equations (via iteration) for steady state 
analyses in which a rub has occurred 

SPEED Processes the rotor speed and rate input (transient and steady state 
analyses) 

SQUEEZ Calculates the forces that the current type 6 physical connecting 

element (squeeze film damper element) exerts on the modal subsystems or 
ground for transient analyses 

STIFFE Computes the engine support-links element stiffness matrix (transient 
and steady state analyses) 

STIFFT Computes the truss (link) stiffness matrix for the trusses (links) that 
are part of the engine support-links element 

SUBSYS Processes the input data for the modal subsystems (transient and steady 
state analyses) 

TILOOP Subroutine for the time integration loop, which finds the solution for 
each time step. This subroutine is the control loop for the solution 
phase for transient analyses. 

TMAT1 Finds the transformation matrix and associated variables for the 
current type 1 physical connecting element (general spring-damper 
element) (steady state analyses only) 

TMAT2 Finds the transformation matrix and associated variables for the 

current type 2 physical connecting element (space link-damper element) 
(steady state analyses only) 

TMAT4 Finds the transformation matrix and associated variables for the 
current type 4 physical connecting element (engine support-links 
element) (steady state analyses only) 

TMAT5 Finds the transformation matrix and associated variables for the 
current type 5 physical connecting element (uncoupled point 
spring-damper element) (steady state analyses only) 

TMATG Finds the transformation matrix and associated variables for the 
current gyroscopic element (steady state analyses only) 

UBAL Processes the input unbalance loads (transient and steady state 
analyses) 



5.0 STEADY STATE ANALYSIS VERIFICATION 


Verification of the steady state option of TETRA 2 is achieved at different 
levels of complexity. This is done partly to check the logic of the problem 
formulation and partly to demonstrate the correctness of the coding. 

The initial stage consisted of checking degenerate cases which have been 
available in the literature for some time. These are usually linear single or 
two degree-of-freedom systems. TETRA 2 results agreed exactly with the 
closed form solutions for these cases (see section 5.1). 

Next, a solution of Duffing’s equation for the nonlinear hardening (or 
softening) spring (rub element) is included (see section 5.2.1). Also, a more 
complex system of a rigid beam with two nonlinear bearings (rub elements) was 
analyzed (see section 5.2.2). For both of these cases, TETRA 2 produced the 
same results as small computer programs which were independently programmed. 

The final test cases use the two component demonstrator model consisting of a 
casing and a rotor. These cases are numerical. Frequency sweep results for the 
demonstator model without any nonlinear springs (rub elements) were compared to 
the results of the NOVAS computer program (a GE program for engine vibrations 
analysis) (see section 5.3.2). 

The demonstrator model including a nonlinear spring (rub element) was also 
analyzed. These results could not be checked against another computer program, 
since we currently do not have any other computer program with this capability. 
However, we did compare the results of the TETRA 2 transient analysis option 
when taken to steady state conditions with the direct solution obtained in the 
steady state option for these cases. This was done for each of the four 
possible rub categories: no rub, continual rub with dead band equal to zero, 

continual rub with dead band not equal to zero, and intermittent rub (see 
section 5.3.3). 

These test cases have verified the logic and the coding of TETRA 2, specifically 
the steady state option. However, the transient analysis option was verified 
also by comparing its results with the original TETRA program. 
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5.1 Degenerate Cases: Closed Form Solutions 


These cases are the classic spring-mass-dashpot systems in one and two degrees 
of freedom. Formulae from references 4 and 5 were programmed independently and 
compared to TETRA 2 results for the simplified models. 

5.1.1 Single Dearee-of-Freedom 

A schematic of the physical model is shown in figure 5-1 and is self 
explanatory. The closed form equations were programmed independently on an IBM 
PC, and the resulting displacements tabulated in figure 5-1. The results using 
TETRA 2 are shown in the displacement versus frequency plot (figure 5-2). As 
can be seen by comparing the tabulated results in figure 5-1 with the plot of 
figure 5-2, exact agreement between TETRA 2 and the closed form equations was 
obtained. 


5.1.2 Two Deoree-of-Freedom 

Figure 5-3 illustrates the other linear degenerate case, where again the 
resulting displacements tabulated in figure 5-3 were obtained by programming 
formulae obtained from references 4 and 5 on an IBM PC. TETRA 2 results for the 
same case are given in figures 5-4 (displacement versus frequency plot for point 
1) and in figure 5-5 (displacement versus frequency plot for point 2). Again, 
as can be seen by comparing the tabulated results of figure 5-3 with the plots 
of figures 5-4 and 5-5, exact agreement between TETRA 2 and the closed form 
equations was obtained. 



C:\TETR A >CF0RM1 


ENTER FO (LB) , K (LB/IN), C ( (LB-SEC) /IN) , AND M (LB) 
1000, 5E5, 87. 8, 500 


ENTER BEGINNING, FINAL, AND STEP W (RPM) 

1000 , 20000, 1000 


FORCING 

DISPLACEMENT 

FREQUENCY 

AMPLITUDE 

« RPM 

Y INCHES 

1000. 

. 00205805 

2000. 

.00225418 

3000. 

. 00267875 

4000. 

. 00363078 

5000. 

.00656400 

6000. 

.01778637 

7000. 

. 00486259 

8000. 

.00241110 

9000. 

.00152751 

10000. 

.00108271 

11000. 

.00081886 

12000. 

. 00064627 

13000. 

. 00052579 

14000. 

. 00043765 

15000. 

. 00037087 

16000. 

.00031886 

17000. 

. 00027744 

18000. 

. 00024384 

19000. 

.00021617 

20000. 

.00019307 


Stop — Program terminated. 



Figure 5-1. Closed Form Equation Results for 1 Degree-of-Freedom Linear 
Test Case 
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Cs \ TETRA >CF0RM2 


ENTER FO (LB), K1 (LB/IN), K2 (LB/IN), C2 < (LB-SEC) /IN) , 
Ml (LB), AND M2 (LB) 

500 , 6E5 , 5E5 , 90 , 1000, 100 

ENTER BEGINNING, FINAL, AND STEP W (RPM) 

1000 , 20000, 1000 


FORCING DISPLACEMENT DISPLACEMENT 

FREQUENCY AMPLITUDE AMPLITUDE 

o RPM Y1 INCHES Y2 INCHES 


1000. 

.00087910 

.00088411 

2000. 

.00105294 

.00107736 

3000. 

.00157399 

.00165844 

4000. 

. 00520835 

.00572521 

5000. 

. 00260335 

.00302909 

6000. 

. 00090976 

.00113924 

7000. 

. 00050937 

.00070014 

8000. 

.00033416 

.00051645 

9000. 

. 00023725 

.00042497 

10000. 

.00017617 

. 00038006 

11000. 

.00013451 

. 00036573 

12000. 

.00010631 

. 00037492 

13000. 

. 00009426 

. 00039070 

14000. 

. 00009849 

. 00035507 

15000. 

. 00009625 

. 00025264 

16000. 

. 00008509 

.00016328 

17000. 

. 00007359 

.00010792 

18000. 

. 00006399 

. 00007466 

19000. 

.00005619 

.00005384 

20000. 

.00004981 

.00004017 

Stop - Program 

terminated. 




Figure 5-3. Closed Form Equation Results for 2 Degree-of-Freedom 
Linear Test Case 
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5.2 Nonlinear Systems: Simple Models 


To assess the iterative procedure for obtaining numerical solutions as well as 
check the logical method for determining the compatibility relations at the 
nonlinear connections, small computer programs were written to analyze the 
following two cases. Iterative subroutines from the IMSL library were initially 
used in these small programs and subsequently employed in TETRA 2. 

Another purpose of writing the small computer programs cases was to obtain an 
independent check of the coding in TETRA 2. The people who programmed these 
small computer programs were not involved in the development or coding of TETRA 
2. When TETRA 2 was run for these test cases, however, the same results as for 
the small computer programs were obtained. This provides confidence that both 
the small computer programs and the very much larger TETRA 2 were coded 
correctly. 

The nonlinearity in the following cases is of the Duffing type, i.e. cubic 
hardening or softening spring with a linear viscous damper. This type of 
nonlinearity is discussed in references 4 and 5. 

5.2.1 Single Dearee-of-Freedom: Nonlinear Spring 

This is the classic Duffing equation with a cubic hardening spring. A sketch 
of the model and the governing equations for this case are shown on the top part 
of figure 5-6. 

A small computer program was written specifically to assess the IMSL library 
subroutine ZSYSTM, which uses Brown’s iteration method. Results of this effort 
were discussed in the monthly progress reports. The most salient conclusion is 
that this subroutine (and also IMSL subroutines ZSCNT and ZSPOW) can yield 
multiple solutions corresponding to the classic "jump" solution of Duffing’s 
equation with large nonlinearity. However, to obtain the solutions in the 
neighborhood (in frequency) of the "jumps", one must run the program twice, once 
with increasing forcing frequency and once with decreasing forcing frequency. 
This is illustrated by the plot of this case shown in figure 5-6. 

This procedure indicates the sensitivity of the iteration routine to initial 
trial solutions. In addition, the use of harmonic balance in the solution of 
Duffing’s equation gives similar results as those obtained by perturbation or 
direct iteration (without harmonic balance). 

TETRA 2 was also originally coded to include IMSL subroutine ZSYSTM. The 
results produced by the small computer program and TETRA 2 were in excellent 
agreement. Soon after this case was run, however, IMSL Incorporated deleted 
subroutine ZSYSTM from their library and replaced it with two new subroutines, 
ZSCNT and ZSPOW. The new iteration subroutines are based on the secant method 
and Powell’s hybrid algorithm (a variation of Newton’s method) respectively. 

Subsequent test cases that we ran confirmed that ZSCNT and ZSPOW both worked 
better than subroutine ZSYSTM, so we also deleted ZSYSTM from TETRA 2 and 
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instead included ZSCNT and ZSPOW. Subroutine ZSYSTM would frequently be unable 
to converge due to a singularity in the Jacobian, which was a great nuisance. 
This problem is not present in the ZSCNT or ZSPOW IMSL subroutines. All three 
IMSL subroutines give the same results when they converge. 
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5.2.2 Planar Rigid Beam with Two Nonlinear Bearings 


This model is the simplest one that checks TETRA 2’s method for obtaining the 
iterating equations, i.e., displacement compatibility at the nonlinear bearings 
(rub elements). A sketch of this model Is shown in figure 5-7. 

The governing equations were obtained two different ways: 

(1) In terms of translation and rotation of the center of gravity 

(2) In terms of two rigid modes: the translation modal distribution 

simulates the cases of pure translation and simple rotation about the 
center of gravity. 

Both formulations gave identical results. The responses that were compared to 
the TETRA 2 output were the physical displacement amplitudes at the bearings 
(rub elements). The first formulation Is summarized here. Following volume 1 
of this report, the generalized coordinate differential equations of motion are: 

Translation. f f (C,4>C t ) !j t Gv.Lj-K.L,)® 

s Ifa)' /**K, *j, * 

Rotation about the center of gravity: - . % 

It * U,L> + U,L. C*l;) 9 tC^L* k.,L}ij iCc^-C, l,)«j 

The symbols are defined in the sketch in figure 5-7. 

These differential equations are transformed to 4 simultaneous nonlinear 
equations by the principle of harmonic balance following the theory of volume 1. 
After some simplification, these become a matrix relating the generalized 
coordinates of the center of gravity on the left hand side to the physical 
displacements and external forces on the right. Thus, 
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Figure 5-7. Planar Rigid Beam with 2 Nonlinear Bearings 
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where: 
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The algebraic matrix equation is first transformed by matrix inversion to obtai 
an explicit expression for the generalized coordinate vector. Letting: 


• Coefficient of 4*,$* 

2 Nonlinear Connecting Force Vector 


■Vt 

^ 2 External Force Vector 

The matrix equation can be written more simply: 


Vr fi * ptu Fi.t fj.-) 

ar dives: ; " w 


Solving for qj gives: 


I , - 


y " 4 



The displacement (cosine and sine) components at the ends may now be expressed 
in terms of and ^ , since: 

A, s a, - djL . 13, sS > " M- 

A t * a, + a t L t ^ + 

Noting that , £*,61 , and $£ are discrete values of^j, and that theCCj are 
expressed as nonlinear functions of A, ,&i, etc., and substituting these into 
the physical displacement expressions gives the iterating equations, thus: 

a , =(e,' l -L,E;;xt i 

A 2 »(%, + + ^i) 

B* •<*£- 

A s -t f 'HO 

The foregoing follows the essential elements of the theory which proceed from 
the formation of the connecting element forces to the formation of a global 
generalized matrix differential equation and to its transformation by harmonic 
balance. It demonstrates also the steps required to obtain the Iterating 
equations from the physical displacement relations at the bearings (rub 
elements). Thus, the application of the theory can be checked, step by step. 

Finally, a small computer program was written to calculate the beam ends’ 
response, and the results agreed. with the TETRA 2 solution. This case also 
revealed potential convergence problems with the IMSL Iteration subroutines. 
Large amplitudes (which amplify the nonlinearities) occasionally led to a large 
number of iterations that would not converge, especially in the regions of 
resonances. However, displacement amplitudes and nonlinear forces more typical 
of actual engine experience avoided the convergence problems. 

For highly nonlinear connecting elements where a "jump” phenomenon is expected, 
manual Inputs of initial Iteration values in TETRA 2 (see input sheet I>2) can 
decrease the number of iterations and overcome non-convergence. 

Figure 5-8 gives some response plots from TETRA 2 for this case. 
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Figure 5-8. Planar Rigid Beam with 2 Nonlinear Bearings 







5.3 Demonstrator Model 

The final set of test cases use the two component demonstrator model consisting 
of a rotor and a casing. This is the exactly the same demonstator model that 
was used back in 1981 to check out the original TETRA program (which was for 
transient analyses only), as documented in reference 1. For more details about 
this model, see reference 1 (especially volume 3, section 7.0). A sketch of the 
model is shown in figure 5-9. 

This demonstrator model was also used in 1983 in modified form (two rotors, 
instead of a rotor and a casing, and the stiffnesses changed) to check out the 
flexible bladed disk option of TETRA (see reference 2). However, the test cases 
made to check out the steady state analysis portion of TETRA conform to the 
original version of the demonstator model (rotor and casing), rather than the 
one that was modified to check out the flexible bladed disk option. 

The earlier test cases for the demonstator model (both with and without flexible 
bladed disk) were repeated, with the same results as before. This verifies that 
the extensive modifications made for TETRA 2 did not mess up the transient 
analysis portion of the program. The transient analysis portion of the program 
was Improved with TETRA 2 by making it more efficient, and by making the new 
features developed for TETRA 2 (such as the improved rub element) applicable to 
both steady state and transient analyses. 

In the following sections, we only concern ourselves with the new test cases 
that were made to verify the steady state analysis portion of TETRA 2. 

Reference 1 (volume 3) and reference 2 give details about the earlier test cases 
that verified the transient analysis portion of the program. To compare with 
the steady state analysis results, some new transient analysis runs were made 
which were taken to state analysis conditions (that is, run for a long enough 
time for the transient analysis part to die out), as reported in section 5.3.3. 

Section 5.3.1 discusses the input and output for a steady state analysis test 
case for the demonstrator model. Then, in section 5.3.2 we compare frequency 
sweep results for the demonstator without any nonlinear springs (rub elements) 
to the results of the NOVAS computer program (GE program for engine vibrations 
analysis). Finally, in section 5.3.3 we examine test cases for the demonstator 
with one rub element, comparing steady state analysis results with transient 
analysis runs which were taken to steady state conditions. As will be seen, 
excellent results were obtained by these test cases, further verifying the new 
steady state analysis capability of TETRA 2. 
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5.3.1 Input and Output 


TETRA 2 input and putput follows for a steady state analysis run of the 
demonstator model. The input listing is a part of the TETRA 2 output and 
immediately follows the title page of the output. This run corresponds to the 
run discussed in section 5. 3. 3. 4 (intermittent rub at 3000 rpm). To produce 
this output, the standard form of the printed output was requested (I0UT - 1 on 
the type A input sheet). 
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THE CURRENT PROGRAM LIMITS FOR THE INPUT VARIABLES ARE: 

MAX NO. OF PHYSICAL POINTS NOT LOCATED ON THE MODAL SUBSYSTEMS= 10 
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INPUT DATA FOR POINTS NOT LOCATED ON THE MODAL SUBSYSTEMS: 
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ROTOR VERTICAL PLANE SUBSYSTEM FOR DEMONSTRATOR MODEL 
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SUBSYSTEM DIRECTIONS: 
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Ai PER INPUT ROUTINE OPTION IROUT = 1, IMSL SUBROUTINE ZSCNT IS USEO 
TO SOLVE THE NONLINEAR EQUATIONS. 
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FORCE AMPLITUDE (POUNDS) 
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PHASE ANGLE FOR THE PHYSICAL VELOCITY 

F'OINI X Y Z THETA-X THETA-Y THETA'Z 

NUMBER DEGREES DEGREES DEGREES DEGREES DEGREES DEGREES 
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THE MAGNITUDES OP THE GYROSCOPIC FORCES ACTInl. ON THE ROTOR(S) AREi 
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5.3.2 Frequency Sweep with No Rub Element Test Cases 

TETRA 2 steady state analysis frequency sweep runs were made for the demonstator 
model with no rub element. Two runs were made, one including gyroscopic 
stiffening at the fan and turbine locations, and one without. For both runs, an 
unbalance of 100 gram- inches was applied at the fan location on the rotor, 
the beginning forcing frequency (rotor speed) was 150 rpm, the step forcing 
frequency 50 rpm, and the final forcing frequency 15000 rpm (a solution 
calculated for each of 298 forcing frequencies). The output for the 298 forcing 
frequencies was written to the output plot file and subsequently plotted, as 
shown in figures 5-10 through 5-13. 

The TETRA 2 results are compared to those produced by the NOVAS computer program 
(GE computer program for engine vibrations analysis), which are shown in figures 
5-14 through 5-17. These figures are the same ones as figures 67 through 
70 of reference 1, volume 3 (the original 1981 TETRA user’s manual). Figure 
5-10 corresponds to figure 5-14 (fan displacement without gyroscopic stiffening 
effects), figure 5-11 corresponds to figure 5-15 (middle of rotor displacement 
without gyroscopic stiffening effects), figure 5-12 corresponds to figure 5-16 
(fan displacement with gyroscopic stiffening), and figure 5-13 corresponds to 
figure 5-17 (middle of rotor displacement with gyroscopic stiffening). As can 
be seen by comparing these figures, excellent agreement between TETRA 2 and the 
NOVAS program was obtained. 
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Figure 5-10. TETRA 2 Steady State Frequency Response for Demonstrator 
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Figure 5-11. TETRA 2 Steady State Frequency Response for Demonstrator 



Figure 5-12. TETRA 2 Steady State Frequency Response for Demonstrator 
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Figure 5-13. TETRA 2 Steady State Frequency Response for Demonstrator 
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Figure 5-14. NOVAS Steady State Frequency Response for Demonstrator 
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Figure 5-15. NOVAS Steady State Frequency Response for Demonstrator 




155 


Figure 5-16. NOVAS Steady State Frequency Response for Demonstrator 
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Figure 5-12. NOVAS Steady State Frequency Response for Demonstrator 



5.3.3 Single Frequency with One Rub Element Test Cases 

Several runs of the demonstrator model with one nonlinear spring (rub element) 
located at the fan were made. One of the prime reasons for these runs was to 
provide a comprehensive checkout of the rub element, including the dead band 
feature of the rub element, which was not exercised in the simpler test cases 
reported in the previous sections. 

For each rub element, TETRA 2 recognizes four possible rub categories: no rub, 

continual rub with dead equal to zero, continual rub with dead band not equal to 
zero, and Intermittent rub. For each forcing frequency, TETRA 2 compares the 
minimum and maximum relative displacement magnitude and the dead band of each 
rub element to determine which of the categories is present for each of the rub 
elements. These four possible rub element categories are illustrated in 
figure 5-18. Knowing which rub category applies is important, because a 
different method is used to find the rub element physical forces for each of the 
four rub element categories. For more details about the four rub categories and 
the methods used for each category, see section 5.3 of volume 1 of this report. 

To fully check out the rub element, demonstator model runs (one steady state run 
and one transient analysis run) were made for each of the four possible rub 
element categories. These runs were made with a rotor speed (forcing frequency) 
of 3000 rpm, and an unbalance of 5000 gm-in applied at the fan location on the 
rotor. The dead band and the spring rates were adjusted to produce a run for 
each of the four rub categories. 

The results when a rub occurs could not be checked against another computer 
program, because we currently do not have any other computer program with this 
capability. However, what we did do is compare the results of the TETRA 2 
transient analysis option when taken to steady state conditions (that is, run 
long enough for the transient analysis part to die out) with the direct steady 
state solution. Since greatly different methods of solution are utilized in the 
transient analysis and the steady state analysis options of TETRA 2, the fact 
that the two methods agree so very well is most encouraging. 
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5.3.3. 1 No Rub 


A steady state analysis test case was run where a rub element located at the 
fan was present but no rub occurred because the dead band which was inputted 
always exceeded the rub element relative displacement magnitude. Two methods 
were used to verify the results of this test case. 

First, since no rub occurred, the rub element was never activated, and the rub 
element physical forces set to zero by the program. Thus, this case should 
produce the same results as for a TETRA run with no rub element, and remember, 
the demonstrator with no rub element was verified by comparing TETRA 2 results 
with that of the NOVAS computer program (see section 5.3.2). A TETRA 2 run with 
the same input except omitting the rub element confirmed that exactly the same 
results are obtained either way. 

The second way that this test case was verified was by comparing the results 
of the steady state analysis run with the results of a transient analysis run 
taken to steady state conditions. The transient analysis run results are shown 
in figures 5-19 through 5-22. The transient analysis run was stopped at .5 
seconds, by which time the transient effect had almost completely died out, as 
can be seen from these plots. 

These results are summarized in table 5-1 below. The steady state analysis 
values in this table are exact numbers from the printout of the run. The 
transient analysis run numbers, however, had to be estimated by reading the peak 
values in the steady state region (.4 to .5 seconds) of figures 5-19 through 
5-22. Although it is hard to read the figures very accurately, it is evident 
that the steady state and transient analysis runs agree quite well. 


Table 5-1. Demonstrator Model Results with No Rub 


Location 


Transient 

Analysis 

Fiaure 

Number 

Rotor at fan station 

.0430 

.043 

5-19 

Case at fan station 

.0287 

.028 

5-20 

Rotor at midpoint station 

.0105 

.0105 

5-21 

Case at midpoint station 

.00869 

.087 

5-22 
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Figure 5-19. Demonstrator Transient Analysis Response Plot 
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Figure 5-20. Demonstrator Transient Analysis Response Plot 
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Figure 5-21. Demonstrator Transient Analysis Response Plot 
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Figure 5-22. Demonstrator Transient Analysis Response Plot 



5. 3. 3. 2 Continual Rub with Dead Band Equal to Zero 

Next, a steady state analysis test case was run for the demonstrator model where 
the rub element, which was located at the fan, had an input dead band equal to 
zero. Since the input dead band was zero, a continual rub of course resulted. 

A transient analysis run taken to steady state conditions was also run to verify 
the results. The transient analysis run results are shown in figures 5-23 
through 5-26. The transient analysis run was stopped at .5 seconds, by which 
time the transient effect had almost completely died out, as can be seen from 
these plots. 

The results are summarized in table 5-2 below. The steady state analysis values 
in this table are exact numbers from the printout of the run. The transient 
analysis run numbers, however, had to be estimated by reading the peak values in 
the steady state region (.4 to .5 seconds) of figures 5-23 through 5-26. 

Although it is hard to read the figures very accurately, again it is evident 
that the steady state and transient analysis runs agree quite well. 


Table 5-2. Demonstrator Model Results for Continual Rub with Dead Band 
Equal to Zero 



Maximum Vertical. 

DisDlacement (inches! 

Transient 

Analysis 

Location 

Steady state 

Analysis Run 

Transient 
Analysis Run 

Fiaure 

Number 

Rotor at fan station 

.0234 

.023 

5-23 

Case at fan station 

.0195 

.019 

5-24 

Rotor at midpoint station 

.00649 

.0066 

5-25 

Case at midpoint station 

.00564 

.0057 

5-26 
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Figure 5-23. Demonstrator Transient Analysis Response Plot 
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Figure 5-24. Demonstrator Transient Analysis Response Plot 


DEMONSTRATOR MODEL UITH GYRO, UITH RUB ELEMENT 
DISPLACW 




167 


Figure 5-25. Demonstrator Transient Analysis Response Plot 
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Figure 5-26. Demonstrator Transient Analysis Response Plot 







5. 3. 3. 3 Continual Rub with De ad Band Not Equal to Zero 

Next, a steady state analysis test case was run for the demonstrator model where 
the rub element, which was located at the fan, had an input dead band not equal 
to zero. The magnitude of the input dead band compared to the minimum and 
maximum relative displacement magnitudes was such that a continual rub resulted 
for this case. A transient analysis run taken to steady state conditions was 
also run to verify the results. The transient analysis run results are shown in 
figures 5-27 through 5-30. The transient analysis run was stopped at .56 
seconds, by which time the transient effect had almost completely died out, as 
can be seen from these plots. 

The results are summarized in table 5-3 below. The steady state analysis values 
in this table are exact numbers from the printout of the run. The transient 
analysis run numbers, however, had to be estimated by reading the peak values in 
the steady state region (.4 to .56 seconds) of figures 5-27 through 5-30. 
Although it is hard to read the figures very accurately, again it is evident 
that the agreement between the steady state and transient analysis runs is 
superb. 


Table 5-3. Demonstrator Model Results for Continual Rub with Dead Band 
Not Equal to Zero 


location 

Maximum Vertical Displacement (inches 1 
Steady state Transient 

Analysis Run Analysis Run 

Transient 

Analysis 

Fiaure 

Number 

Rotor at fan station 

.0303 

.030 

5-27 

Case at fan station 

.0227 

.022 

5-28 

Rotor at midpoint station 

.00789 

.0079 

5-29 

Case at midpoint station 

.00672 

.0067 

5-30 


Also shown are some orbit plots for this case, which are helpful in visualizing 
what is happening. The first orbit plot, figure 5-31, shows the vertical versus 
the horizontal displacement of the points located on the rotor and case at the 
fan station (these points are also the joints of the rub element). This plot is 
for the transient analysis run, and covers through .05 seconds of the run in 
order to show clearly the early transient portion of the run. The displacements 
are 0 at time - 0 (time ■ 0 represents the time of fan blade loss), then 
displacements increase as is shown by the movement out from the center of the 
plot. As can be seen from the plot, the rotor and case points move almost 
parallel to one another, with the rotor displacement larger than the case 
displacement. 
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FIGURE 5-27. DEMONSTRATOR TRANSIENT ANALYSIS RESPONSE PLOT 


DEMONSTRATOR MODEL UITH RUB ELEMENT, TRANSIENT ANALYSIS 



FIGURE 5-28. DEMONSTRATOR TRANSIENT ANALYSIS RESPONSE PLOT 




The next orbit plot, figure 5-32, shows only the point on the rotor at the fan 
station, but for the full 0 through .56 seconds of the transient analysis run. 
Likewise, figure 5-33 shows only the point on the case at the fan station for 
the full 0 through .56 seconds of the transient analysis run. Displacements 
start at 0 at time - 0 (the center of the plot), then there is an overshoot 
where for just a revolution or two the displacements are appreciable higher than 
the rest of the run, after which the motion settles into a steady state phase at 
a constant radius. This is the same pattern exhibited by figures 5-27 through 
5-30. 

Next, figure 5-34, like figure 5-32, shows only the point on the rotor at the 
fan station. Likewise, figure 5-35, like figure 5-33, shows only the 
point on the case at the fan station. However, figures 5-34 and 5-35 were 
generated from the steady state analysis run, whereas figures 5-32 and 5-33 were 
generated from the transient analysis run. Since figures 5-34 and 5-35 plot 
steady state motion only, the transient effects of displacement build-up and 
overshoot are gone, and what is left is a circular orbital path. By overlaying 
figure 5-34 on top of figure 5-32, and figure 5-35 on top of figure 5-33, we can 
see that the steady state analysis results agree very well with the transient 
analysis results after the transient effect has died out. 

Figure 5-36 shows a different type of orbit plot. This plot is for the 
transient analysis run for the full 0 through .56 seconds. However, this plot 
shows the relative displacement components for the rub element on the x and y 
axes, whereas the earlier orbit plots showed the actual physical displacement 
components of the two joints of the rub element on the x and y axes. (The 
relative displacement for the rub element is defined as the displacement of 
joint I of the rub element minus that of joint J). One thing that is' neat about 
this type of orbit plot is that we also show the clearance circle on this plot. 
The radius of the clearance circle is the dead band which was inputted, in this 
case .005 inches. As can be seen, the dead band is exceeded (it rubs) for all 
except the very early portion of the transient analysis run. 

Finally, figure 5-37 gives the same type of orbit plot as figure 5-36, except 
that figure 5-37 was generated from the steady state analysis run, whereas 
figure 5-36 was generated from the transient analysis run. Again, by overlaying 
figure 5-37 on top of figure 5-36, we can see that the steady state analysis 
results agree very well with the transient analysis results after the transient 
effect has died out. 

Note that this case produced circular orbital paths (at least after the 
transient effect died out). The reason that circular orbital paths are produced 
is that the input stiffnesses of the physical connecting elements are the same 
in the vertical and horizontal directions, so neither direction has a tendency 
to displace more than the other under the influence of the unbalance load. 
Incidentally, the cases discussed in sections 5.3.3. 1 and 5. 3. 3. 2 also had input 
stiffnesses that were the same in the vertical and horizontal directions, and 
those cases would also produce circular orbital plots similar to the ones shown 
in this section. 
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FIGURE 5-32. DEMONSTRATOR TRANSIENT ANALYSIS ORBIT PLOT 
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FIGURE 5-34. DEMONSTRATOR STEADY STATE ANALYSIS ORBIT PLOT 
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FIGURE 5-35. DEMONSTRATOR STEADY STATE ANALYSIS ORBIT PLOT 
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5. 3. 3. 4 Intermittent Rub 

Finally, a steady state analysis test case was run where an intermittent rub 
resulted. An intermittent rub occurs less frequently than no rub or a continual 
rub, and the input had to be adjusted carefully to make it happen. The previous 
demonstrator test cases had the same spring constants in the vertical and 
horizontal directions. However, this results in circular orbital paths (see 
figures 5-31 through 5-37). Circular orbital paths will not cause an 
intermittent rub, because the clearance circle is also circular, and a circular 
orbital path will either lie within the clearance circle (no rub) or outside the 
clearance circle (continual rub). 

To cause an intermittent rub, the mount springs for the casing (springs 1 and 2 
in figure 5-9) were made 5 times more stiff in the vertical direction than in 
the horizontal direction. Stiffer casing springs in the vertical direction than 
the horizontal direction might result, for instance, from an engine attached to 
a wing via a pylon. This discrepancy in spring rates in the two directions 
forced the orbital path to take an elliptical shape. The input dead band was 
then adjusted so that the orbital ellipse intersected the clearance circle at 
four points, causing an intermittent rub, as illustrated in figure 5-18. 

As for the previous test cases for the demonstrator model with one rub element, 
a transient analysis run taken to steady state conditions was made for 
comparison. The transient analysis results are shown in figures 5-38 through 
5-45. The transient analysis run was stopped at .56 seconds, by which time the 
transient analysis effect had almost died out, as can be seen from these plots. 
Plots for the horizontal direction displacements are included a well as those 
for the vertical. Only the vertical direction displacement plots were included 
for the test cases described in sections 5.3.3. 1, 5. 3. 3. 2, and 5. 3. 3. 3, because 
the vertical and horizontal displacement plots were essentially the same for 
those test cases due to the circular orbital paths.. 

The results are summarized in table 5-4 on page 191. For the test cases 
described in sections 5.3.3. 1, 5. 3. 3. 2, and 5. 3. 3. 3, the maximum horizontal 
displacements and maximum vertical displacements were the same since the orbital 
paths were circular. Therefore, only the maximum vertical displacement results 
were tabulated for those test cases. However, for the intermittent rub test 
case the orbital path is elliptical, and the maximum horizontal displacements 
are much larger than the maximum vertical displacements. Therefore, both 
vertical and the horizontal displacement results are summarized in table 5-4. 

The steady state analysis values In this table are exact numbers from the 
printout of the run. The transient analysis run numbers, however, had to be 
estimated by reading the peak values in the steady state region (.4 to .56 
seconds) of figures 5-38 through 5-45. 

This test case provides the most difficult challenge for the steady state 
analysis to meet. As is seen from table 5-4, agreement between the steady state 
and transient analysis runs is very good for most point and directions. However, 
agreement is not as good as for the test cases of sections 5. 3. 3.1, 5. 3. 3. 2, and 
5. 3. 3. 3, and agreement is poor for one point and direction (midpoint of rotor in 
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Figure 5-38. Demonstrator Transient Analysis Response Plot 




Figure 5-39. Demonstrator Transient Analysis Response Plot 
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Figure 5-41. Demonstrator Transient Analysis Response Plot 
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Figure 5-42. Demonstrator Transient Analysis Response Plot 
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Figure 5-43. Demonstrator Transient Analysis Response Plot 
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the vertical direction). 


Table 5-4. Demonstrator Model Results with Intermittent Rub 


Location 

Direction 

Maximum Displacement (inches) 
Steadv State Transient 
Analysis Run Analysis Run 

Transient 

Analysis 

Fiaure 

Number 

Rotor at fan station 

Horizontal 

.03905 

.039 

5-38 

Rotor at fan station 

Vertical 

.00479 

.0045 

5-39 

Case at fan station 

Horizontal 

.02961 

.030 

5-40 

Case at fan station 

Vertical 

.00149 

.0017 

5-41 

Rotor at midpoint station 

Horizontal 

.00698 

.0067 

5-42 

Rotor at midpoint station 

Vertical 

.00148 

.0030 

5-43 

Case at midpoint station 

Horizontal 

.00596 

.0060 

5-44 

Case at midpoint station 

Vertical 

.00076 

.0008 

5-45 


The reason that this test case does not agree as well as the previous test cases 
has to do with the interaction between the rotor and case resulting from the rub 
element. For a transient analysis run, the rub element works as follows. For a 
good portion of the orbit, the clearance circle is not exceeded, and the rub 
element provides no force. However, twice during each orbit the clearance 
circle is exceeded and a rub element force is exerted on the rotor and case. 

This on/off rub force has a complex effect on the vibration of the rotor and the 
case. Rather than just first harmonic vibration, involving the terms coswt and 
sin<ot, higher order terms involving cos2ut, sin2cot, cos3<ot, sin3a>t, etc. may be 
produced . 

Looking at the transient analysis response plots for this case, figures 5-38, 
5-39, 5-40, and 5-44 appear to involve only simple first harmonic vibration, 
with a frequency of vibration oi equalling the forcing frequency, which in this 
case is the rotor speed (3000 rpm). However, in figures 5-41, 5-42, 5-43, and 
5-45, two frequencies of vibration are observed - the forcing frequency of 3000 
rpm and also three times that (9000 rpm). So, in addition to first harmonic 
vibration, involving the terms cosat and sinwt, we have third harmonic 
vibration, involving the terms cos3<ot and sin3wt. 

Steady state analyses are different. For steady state analyses, the technique 
of harmonic averaging is used to find the rub element forces, as explained in 
section 5.2 of volume 1. Rather than apply the rub force in an on/off manner, 



for steady state analyses the rub force is smeared over the entire orbit. All 
higher order terms are neglected as part of harmonic averaging, leaving only the 
first harmonic terms. For this reason, for intermittent rubs the transient 
analysis taken to steady state conditions provides a more accurate 
representation of the true characteristics than does the steady analysis. 

This is also the reason for the poor agreement between the steady state and 
transient analyses for the rotor at the midpoint station in the vertical 
direction. The transient analysis response plot for the rotor at the midpoint 
station in the vertical direction (figure 5-43) shows that the 9000 rpm mode 
(the third harmonic mode) is very active for this particular point and 
direction. Herein lies the rub (pardon the pun), since this mode is completely . 
neglected in the steady state analysis. This was the only one of the eight 
transient analysis response plots where this mode was so very strong compared to 
the first harmonic mode, which explains why this point and direction was the 
only one to have poor agreement. Note that this point and direction are of very 
little importance, since the vibration levels are small compared to the 
vibration in the horizontal direction. 

Orbit plots for this case are included in figures 5-46 through 5-52. These 
orbit plots are similar to the ones in section 5. 3. 3. 3, except that the orbital 
paths are elliptical in shape rather than circular. Section 5. 3. 3. 3 discusses 
the characteristics of the orbit plots in detail. 

The steady state analysis orbital paths for this case agree very well with the 
transient analysis orbital paths after transient effects have died out, as can 
be seen by comparing figure 5-49 with figure 5-47, figure 5-50 with figure 5-48, 
and figure 5-52 with figure 5-51. One difference, however, is that the steady 
state analysis orbital paths are true ellipses, whereas the transient analysis 
orbital paths after the transient effect has died out are somewhat distorted 
ellipses. This distortion results from the complex interaction between the 
rotor and case transmitted through the rub element, with the rub element force 
being activated when the orbit exceeds the clearance circle, and deactivated 
when inside the clearance circle. The transient analysis orbit plots, of 
course, are more representative of the true orbital paths, because the steady 
state analysis is an idealization which does not account for the full complexity 
of the rotor and case interaction. 

In summary, excellent results were obtained in verifying the steady state 
analysis option of the TETRA 2 computer program. Exact agreement between TETRA 
2 and closed form equations was obtained for simple linear one and two 
degree-of- freedom cases. Also, for simple nonlinear one and two 
degree- of- freedom test cases (with one and two rub elements, respectively), 
excellent agreement was obtained between TETRA 2 and small computer programs 
which were written independently. Using the demonstrator model, very good 
agreement was obtained between TETRA 2 and the NOVAS computer program for a 
frequency sweep with no rub elements. Finally, for the demonstrator model with 
one rub element, very good agreement between the steady state and transient 
analysis options of TETRA 2 was obtained for each of four possible rub 
categories. 
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Figure 5-46. Demonstrator Transient Analysis Orbit Plot 









Figure 5-47. Demonstrator Transient Analysis Orbit Plot 
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Figure 5-48. Demonstrator Transient Analysis Orbit Plot 
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Figure 5-49. Demonstrator Steady State Analysis Orbit Plot 



Figure 5-50. Demonstrator Steady State Analysis Orbit Plot 




DELTA X HORIZONTAL, INCHES 







6.0 REFERENCES 


1. Gallardo, V. et al, ’Blade Loss Transient Analysis’, Task II and Task III, 
NASA CR- 165373, Volumes II and III, June 1981. 

2. Gallardo, V. and Black, G. et al, ’Blade Loss Transient Dynamic Analysis 
with Flexible Bladed Disk’, NASA CR- 168176, April 1983. 

3. Ghaby, R.A., ’The Transient/Nonlinear Vibration of Gas Turbine Jet Engines 
With Squeeze Film Dampers Due To Blade Loss’, Master’s Thesis, Department 
of Mechanical and Aerospace Engineering, Case Western Reserve University, 
Cleveland, Ohio, May 1984. 

4. Dinca, F. and Teodosiu, C., ’Nonlinear and Random Vibrations’, Academic 
Press, N.Y., N.Y., 1973. 

5. Minorsky, N., ’Introduction to Nonlinear Mechanics’, J.W. Edwards, Ann 
Arbor, MI., 1946. 


200 



